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Preface 



It is known that the balance between nonlinearity and dispersion may result 
in an appearance of localized long bell-slraped strain waves of permanent 
form (solitary waves or solitons) which may propagate and transfer energy 
over the long distance along. Starting with the first documented water 
surface solitary wave observation, made by J. Scott Russell yet in 1834, 
solitons in fluids were observed and generated many times. It was the 
most surprising fact, however, that despite of almost similar description of 
stresses in fluids and solids, bulk longitudinal strain solitary waves have not 
been observed in nonlinearly elastic wave guides. 

One of reason of the lack of the results on nonlinear wave in solids is 
that the complete description of a three-dimensional (3-D) nonlinear con- 
tinuum is a difficult problem. That is why initial 3-D problems are usually 
reduced to the one-dimensional (1-D) form in order to clarify the simplest 
but qualitatively new analytical solutions. Certainly the cylindrical elastic 
rod seems to be a suitable real-life 1-D wave guide. Recently, the theory has 
been developed to account for long longitudinal strain solitary waves prop- 
agating in a free lateral surface elastic rod with permanent cross section. 
The nonlinearity, caused by both the finite stress values and elastic mate- 
rial properties, and the dispersion resulting from the finite transverse size 
of the rod, when in balance allow the propagation of the bulk strain soli- 
tary waves. Motivated by analytical theoretical predictions, there has been 
successful experimental generation of strain solitary waves in a polystyrene 
free lateral surface rod using the holographic interferometry. Hence it was 
proven that bulk long localized nonlinear strain waves of permanent form 
really exist. 

However, presence of a dissipation (accumulation) destroys the balance 
between nonlinearity and dispersion, and nonlinear strain wave in the rod 
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viii Book Title 

may attenuate or amplify. One possibility occurs when the radius of the 
rod varies. Dissipative (active) effects may be caused by internal features 
of the elastic material, hence, an irreversible part should be included into 
the stress tensor in addition to the reversible one depending only upon the 
density of the Helmholtz energy. Dissipation (accumulation) may also come 
in an elastic wave guide through phenomena occurring at its lateral surface. 
Among the volume sources of dissipation/accumulation one can mention a 
microstructure, point moving defects, thermal effects. When dissipation 
and activation act together there may be another balance resulted in a for- 
mation of a bell-shaped wave with the amplitude and the velocity prescribed 
by the condition of the dissipative/active balance. Hence the wave is se- 
lected. Note that there exist another kind of nonlinear wave of a permanent 
shape sustained either by a balance between nonlinearity and dissipation or 
by a balance between nonlinearity, dispersion and dissipation. This wave 
has the form of a shock and is often called kink-shaped wave or simply kink. 

The amplification of the waves (i.e., growth of the amplitude) may cause 
the appearance of plasticity zones or microcracks in a wave guide. This is 
of importance for an assessment of durability of elastic materials and struc- 
tures, methods of non-destructive testing, determination of the physical 
properties of elastic materials, particularly, polymeric solids, and ceramics. 
Bulk waves provide better suited detection requirements than surface strain 
waves in setting up a valuable non-destructive test for pipelines. 

Inclusion of dissipation (accumulation) yields nonlinear dispersive- dis- 
sipative governing equations that are nonintegrable as a rule. Hence, only 
particular, usually travelling wave, solutions may be obtained analytically. 
Certainly, these solutions require specific initial conditions. Moreover, they 
usually have no free parameters, and additional relationships between the 
equation coefficients are required for the existence of the solutions. That is 
why the obtaining of the exact solutions is often considered as useless by 
many authors preferring to apply numerical methods only. 

I would like to achieve two tasks in this monograph. First, it is planned 
to provide the sequential analytical consideration of the strain waves am- 
plification/attenuation and selection in solids, mainly in an elastic rod. It 
may be of interest for those working in the field of solid mechanics. Another 
task is to demonstrate the use of even particular analytical solutions for the 
description of unsteady nonlinear wave processes. It may attract the atten- 
tion of the specialists in various fields since the structure of the governing 
equations is rather universal. The content is essentially based on the author 
previous research. However, many works were done in a collaboration. The 
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Chapter 1 

Basic concepts 



This chapter is focused on some features of nonlinear waves to be used 
further in the book. Linear waves are accounted for the linear equations, 
they have infinitesimal amplitude. Nonlinear waves are described by non- 
linear equations. In contrast to the linear waves, an amplitude, a veloc- 
ity and a wave number of the nonlinear waves are connected to one an- 
other. More general information about nonlinear waves may be found in 
numerous special books, like Ablowitz and Segur (1981); Bhatnagar (1979); 
Calogero and Degasperis (1982); Newell (1985); Sachdev (1987); Whitham 
(1974) etc. 

The governing equations for the nonlinear strain waves to be considered 
are nonitegrable by the inverse scattering transform method, and only par- 
ticular exact solutions may be obtained. Of special interest are the single 
travelling wave solutions that keep their shapes on propagation. This is a 
result of the balances between various factors affecting the wave behaviour. 
There are two main types of the nonlinear travelling solitary waves which 
could propagate keeping its shape, bell-shaped and kink-slraped solitary 
waves. The bell-slraped solitary wave usually appears as a result of a bal- 
ance between nonlinearity and dispersion. The kink-slraped wave may be 
sustained by different balances, one possibility occurs when nonlinearity is 
balanced by dissipation (or accumulation) , another case corresponds to the 
simultaneous balance between dispersion, nonlinearity and dissipation (or 
accumulation). The single travelling wave solution requires specific initial 
conditions. However, one can show that these solutions account for the 
final quasistationary part of an arbitrary initial pulse evolution. This un- 
steady process may be described analytically for the integrable equations 
or numerically for others. We illustrate all mentioned above further in this 
Chapter. 



1 
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Amplification of Nonlinear Strain Waves in Solids 



1.1 Single nonlinear waves of permanent shape 

1.1.1 Monotonic bell- shaped solitary waves 

The simplest celebrated model equation containing nonlinear and dispersive 
terms is the well-known Korteweg- de Vries (KdV) equation Korteweg and 
de Vries (1895), 

u t + 2bu u x + du xxx = 0, (1.1) 

whose exact one-parameter single solitary wave solution is 

u = 6yk 2 cosh -2 k(x — Mk 2 t). (1.2) 

b 

The wave amplitude A = 6dk 2 /b and the velocity V = Adk 2 depend 
upon the wave number k which is a free parameter. One can call the solution 
(1.2) travelling solitary wave one since it depends upon the phase variable 
9 = x—Vt only, and monotonic solitary wave since it decays monotonically 
when |0| — > oo. Typical shape of the wave is shown in Fig. 1.1 where one 
can see also that the wave is symmetric with respect to its maximum. 




Fig. 1.1 Monotonic solitary wave (solid line) and its first derivative (dashed line) 

Sometimes there is a need for the inclusion of higher- order deriva- 
tive (dispersion) or nonlinear terms into Eq.(l.l). A particular case arises 
for water waves when surface tension suppresses coefficient d Hunter and 
Sclreurle (1988) and fifth-order dispersion u$ x is added in Eq.(l.l). 
Also higher- order derivative terms model weak nonlocality Engelbrecht 
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and Braun (1998), provide an improvement of bad dispersive properties 
Christov et. al (1996); Maugin and Muschik (1994), account for a contin- 
uum limit of discrete models with far neighbour interactions Koseviclr and 
Savotchenko (1999), to say nothing of dissipative (active) generalizations. 
An example of the inclusion of higher- order nonlinearity is the Sawada- 
Kotera equation Sawada and Kotera (1974). 

Hence, the following nonlinear equation may be considered: 

u t + 2 buu x + 3cu 2 u x + ru u xxx + su x u xx + du 3x + f u 5x = 0, (1.3) 

We get from Eq.(1.3) a fifth-order (in space derivatives) KdV equation 
Hunter and Sclreurle (1988) when c = r = s = 0. This equation was studied 
in many papers, see, e.g., Karpman (1993); Karpman (1998); Karpman and 
Vanden-Broeck (1995); Kawahara (1972); Kawahara and Takaoka (1988); 
Benilov et. al (1993); Grimshaw et. al (1994). When, in addition d = 0, 
the resulting equation models the LC ladder electrical transmittion lines. 
Its solutions were obtained in Nagashima and Kuwahara (1981); Kano and 
Nakayama (1981). A special integrable case corresponds to the Sawada- 
Kotera equation with b = d = 0, c = — r = — s/2 = 10, / = 1 Sawada 
and Kotera (1974). Its solitary wave solutions may be found in Parkes and 
Duffy (1996), see also references therein. 

Equation (1.3) is obviously nonintegrable by the Inverse Scattering 
Transform method, and only particular exact solutions may be obtained. 
Let us consider an exact solution vanishing at infinity. In case of the fifth- 
order KdV equation it has the form Kano and Nakayama (1981): 

u — ^ k 2 cosh~ 4 k{x — Vt), (1.4) 

loo 

with k 2 = — d/(52f). Hence the width of the wave is prescribed by the 
dispersion coefficients d and f which should be of opposite sign. For the 
wave velocity we have V = 144dfc 2 /13 = — 36d 2 /(169/). Hence, simultane- 
ous triggering of the signs of b, d and / results in changing only the wave 
propagation direction. 

In the general case the exact solitary- wave solution has a form similar 
to the KdV soliton (1.2), 



u = A cosh 2 k( x — Vt). 



(1.5) 



Important particular cases are: 
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Amplification of Nonlinear Strain Waves in Solids 



(i) In presence of only cubic nonlinear term, r — s = 0, we get the solution 
with fixed parameters, 



30/ 



A = 2\ J -k 2 , V = Ak z {2b\ — — 5d— 116/ k 2 ), k 2 = 



30/ 



b^/— 30 f/c — 3d 
60/ ' 



(1.6) 

In this case the existence of solution vanishing at infinity is provided by 
the linear fifth order term f u 5x . Indeed at / = 0 we get from Eq.(1.3) the 
Gardner equation whose solution is Grimshaw et. al (1999): 



A\ 

cosh m( x —Vt)+Bi_' 



(1.7) 



where 



. 3\f 7 2drn 2 \[2b , 2 , . 

Ai = , Bi = , V = dm 2 . (1.8) 

V 2b 2 + 9 cdm 2 \j2b 2 + 9 cdm 2 



An important feature of the solution is the existence of the finite limiting 
amplitude when B is large Slyunyaev and Pelinovsky (1999). However, at 
nonzero / a substitution of Eq.(1.7) into Eq.(1.3) yields Bi = 1, cosh m( x— 
Vt) + 1 = 2cosh 2 m(a; — Vt)/ 2, and we get the solution (1.5), (1.6) with 
m = 2k. 

(ii) When only quadratic higher- order term r uu xxx is taken into account, 
c = s = 0, the fixed parameters of the solution (1.5) are 



A _ 30/fc 2 T/ _ 2(50 6 2 / 2 + 5 bdfr - 3 d 2 r 2 ) ; 2 _ 56/ - dr n ^ 
r ’ 25 f r 2 ’ 10 fr ' 1 j 

Note that the solution exists at d = 0. 

(iii) In case c = r = 0 the parameter k is free but an additional restriction 
on the equation coefficients holds, 

A= 6 °^,E = 4fc 2 (d + 4/A: 2 ),106f = ds. (1.10) 

s 

One can see that s may be excluded from the amplitude expression using 
the third formula from (1.10). Then the amplitude coincides with that of 
the KdV soliton (1.2). We also see that the wave velocity consists of two 
parts, V\ = 4 dk 2 and E 2 = 16/A: 4 , the first of which being exactly the KdV 
soliton velocity. Let us rewrite the ODE reduction of the equation (1.3) in 
the form ( ' = d/dd , 9 = x — Vt): 



bu 2 + d u" — V\u + ^ u' 2 + f u"" — V 2 u = 0. 




February 11, 2004 20:10 



WSPC/Book Trim Size for 9in x 6in 



ws-book9x6 



Basic concepts 



5 



One can check that the solution (1.5), (1.10) satisfies separately 



bu 2 + d u" — V\u = 0, (1.11) 

S -u ,2 + fu""-V 2 u = 0, 



where the first of these equations is the ODE reduction of the KdV equation 
(1.1) having a one-parameter solitary- wave solution. 

(iv) Higher order nonlinear terms may support the existence of solitary- 
wave solutions even in absence of the linear dispersion terms. Higher- order 
nonlinear terms provide bounded localized solutions at d = 0 in contrast to 
the case c = r = s = d= 0 Kano and Nakayama (1981). The parameters 
of the solution (1.5) are 



A = 



120 fk 4 
b + 2 k 2 {r + s) 



, V = 16 fk 4 



while k satisfies the equation 



4[30c/ — r(r + s)]/c 4 + 2 bsk 2 + b 2 = 0. 



Thanks to the higher order terms the solitary- wave solution may exist 
even in the absence of the KdV’s nonlinear term, 6 = 0, provided that the 
restriction 30c/ — r(r + s) = 0 is satisfied. Then k may be a free parameter. 

When / = 0 we have 

^ 3 cd — 2b(2r + s) 2d[3cd — 2b(2r + s)] 3cd — 2b(2r + s ) 

c(r + s) ’ (r + s) (2 r + s) ’ 2(r + s) (2 r + s) ' 

There is no exact solution vanishing at infinity in the case d = / = 0. 
Instead the solution in the form of a solitary wave on an ’’pedestal” may 
be obtained as 



u = Hcosh 2 k(x — Vt) + B, (1.12) 

with 

2k 2 {2r + s) (2r + s)[b + 2k 2 (r + s)] 

= , Lj = , 

c 3c(r + s) 

s(2r + s)[4fc 4 (r + s) 2 — b 2 ] 

3 c(r + s) 2 

Even equations with dissipation may possess bell-shaped solitary wave 
solution. In particular, it was recently found Garazo and Velarde (1991); 
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Rednikov et. al (1995) that appropriately heating a shallow horizontal liq- 
uid layer long free surface waves u(x, t) can be excited whose evolution is 
governed by a dissipation-modified Korteweg- de Vries (DMKdV) equation 



lit “t“ 2oq U 'U x ~ t~ OL2 U X x T ^xxx~k ^4 ^xxxx~k ^5 — 0- (1.13) 

The coefficients in Eq.(1.13) depend upon parameters characterizing the 
liquid (Prandtl number etc.), temperature gradient across the layer, and its 
depth. The exact travelling bell-shaped solitary wave solution have been 
obtained in the form (1.12) Lou et. al (1991); Porubov (1993) with 

A = 12 a 4 k 2 / a$, B = — («2 +4 a 4 k 2 )/ a$, 



V = — 2oq a .2 / cx 5, «3 = 2oi\ a 4 / 015 . 



(1.14) 



The meaning of the last expression in (1.14) is similar to that in case c = 
r = 0 for Eq.(1.3). Indeed, when the relationships for V and a 3 hold, the 
ODE reduction Eq.(1.13) may be rewritten as 



/ ~ & 2g;i Ot2 x 0^4 ^5 / 2\ P) / -1 I r\ 

(a 2 wx H ^ ) Ug+— u ss9 -\ —(u)g =0, (1.15) 



dd 



a 5 



Q!2 



2 «2 



The restrictions on the equation coefficients do not necessary provide 
an evidence of the KdV ODE reduction like Eqs.(l.ll), (1.15). Particular 
case corresponds to the Kawahara equation (Eq.(1.13) with a 5 = 0) whose 
exact solution is Kudryashov (1988): 

1 p- 3 

u = 2 coslr~ 2 ( 6>)(1 — tanlr( ^ 6 >)), (1-16) 

128«i a 4 8 a 4 8 a 4 



where V = 5a|, a^= a|/(16 a 4 ). 

All solutions (1.4), (1.5), (1.7) account for monotonic and symmetric 
solitary waves. Despite difference in their functional form they have one and 
the same shape shown in Fig. 1.1. In contrast to them exact solitary wave 
solution (1.16) of the Kawahara equation is monotonic but asymmetric, see 
Fig. 1.2. 



1 . 1.2 Oscillatory bell-shaped solitary waves 

The solitary wave does not decay necessarily in a monotonic manner. Thus 
Kawahara (1972) studied decay at infinity of the wave solution of the fifth- 
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u 



b 




Fig. 1.2 Symmetric (solid line) vs asymmetric monotonic solitary wave (dashed line), 
a) profiles; b)their first derivatives. 



order KdV equation using linearized equation analysis phase analysis. It 
was found that the wave decays monotonically or oscillatory depending 
upon the parameter e, which is proportional (in our notations) to d and 
is inverse proportional by the product of / and the wave velocity. The 
same technique has been used in Karpman (2001) when nonlinear term in 
the fifth-order KdV equation is of the form u p u x . The oscillatory solitary- 
wave solution is shown in Fig. 1.3. The profile of the first derivative with 
respect to the phase variable reveals its symmetric nature. 

Eq.(1.3) does not possess exact oscillatory solitary wave solution. How- 
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u 




Fig. 1.3 Oscillatory solitary wave (solid line) and its first derivative (dashed line) 

ever, it may be described asymptotically. Certainly the fifth-order KdV 
equation is often considered as a perturbed KdV equation. First the asymp- 
totic solution is obtained which consists of the KdV solitary wave solution 
and a small perturbation that oscillates but does not vanish at infinity 
or a non-local solution Benilov et. al (1993); Hunter and Sclreurle (1988); 
Karpman (1993); Karpman (1998). Let us consider the fifth- order deriva- 
tive term and higher- order nonlinear terms as small perturbations assum- 
ing / = SF , r = SR, s = SS, c = SC, S « 1. The asymptotic solution 
u = u{9), 9 = x — Vt, is sought in the form 

u = u 0 (9) + 5 ui(9) + ..., (1-17) 

with Ui — y 0 at \9\ — > oo. In the leading order we get the KdV equation 
(1.1) for the function uo whose travelling solitary wave solution is (1.2). In 
the next order an inhomogeneous linear equation results for u\, 

2b(u 0 ui)g + dui t eee - 4 dk 2 = -F u 0t ggegg - 3C itpMo.e - Ru 0 uo,ggg- 



S UQfiUofig. 



(1.18) 



Its solution vanishing at infinity is 



Ml = ( 



3C 

~4b 



5 bF 
2cP 



2 R + S 




4 R 
~b 



14 F 
~cT 



9 Cd 



)k' 2 M 0 



2 k 2 F 
cL 



9u 0 ,g. 

(1.19) 



b 2 
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The shape of the solution u = u 0 (9) + S Ui(6) depends upon the values 
of the coefficients of Eq.(1.3). It may account for an oscillatory solitary 
wave solution first obtained numerically in Kawalrara (1972), see Fig. 1.3. 
In the case of the fifth-order KdV equation, this profile exists at positive d 
and /, what corresponds to the Case IV in Kawalrara (1972). 

Finally, there exist nonlinear equations having exact travelling wave 
solutions in the form of an oscillatory solitary wave. In particular, an 
equation 

Ut + au x + 2buu x + 3cu 2 u x + du^ x + f u$ x + guj x = 0, (1.20) 

possesses the exact solution 

u = . r^Lk 3 cosh -1 k(x — Vt) 

V c 

where 6 = 0, k 2 = 17//(581g), V = a + 102825fc 6 g/289. The equation 
coefficients should be connected by dg = 37.405 f 2 . 

Oscillatory solitary waves of permanent shape arise also in dissipative 
problems, however, usually they are asymmetric and may be found only 
numerically. 



( 288 \ 
24coslffi 2 k(x — Vt) — — J , 



1 . 1.3 Kink-shaped waves 

The celebrated Burgers’ equation Burgers (1948) is the simplest equation 
that models the balance between nonlinearity and dissipation, 

u t + 26 uu x + gu xx = 0, (1-21) 

In particular, it possesses the exact travelling solution of permanent form, 



u = A m tanh(m 9) + B, 



(1.22) 



where 

g V 

A= B = — , to- free. 
If the boundary conditions are 



— oo, 



(1.23) 
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then 



m = 



(hi - h 2 )b 

2 9 



V = 6(/ii + h 2 ). 




Fig. 1.4 Burgers’ kink-shaped wave 

The shape of the solution (1.22), called kink, is shown in Fig. 1.4. Kinks 
may arise also due to the balance between nonlinearity, dispersion and dissi- 
pation like in the case of the Korteweg-de Vries-Burgers equation (KdVB), 

u t + 2buu x + gu xx + du xxx = 0, (1-24) 

whose exact solution was obtained independently by many authors Vlieg- 
Hultsman and Halford (1991) 

u = A tanh(m 0)sech 2 (m 9) + 2A tanh (m9) + C, (1.25) 

with 

A _ 6 g 2 r _Y_ „ _ 9 

50 Vd' 2 b’ 10d 

It follows from the boundary conditions (1.23) that 
h + — h- = 2 B, V = b(h+ + h _), 
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and the solution exists under 



h 2 -h 2 - —L 
h+ - ~ 25 bd' 



Equations without dissipative terms may also possess the kink-shaped 
solutions. Thus the modified Korteweg-de Vries equation (MKdV), 



Ut “t - 3 cu u x T du xxx — 0 , (1-26) 

has an exact solution 



u= y -to tanh(m0), (1-27) 

where V = —dm 2 . Note that the MKdV equation has both the bell-slraped 
and the kink-slraped solutions. Dissipative equation may possess the same 
property. In particular, the DMKdV equation (1.13) has exact kink-slraped 
solution Lou et. al (1991); Porubov (1993) 



u = C + — ^ — —to cosh 2 (\/~ 3 md)+ 

ot 5 



12(03 S 5 -2 Ql S 4 ) y=3^ tanh( y=3^ 0)) 
5 o 5 



(1.28) 



where 



12 04 02 201(03 05 — 2oi 04) (03 05 — 2oi 04) 2 

G = — — to — — -I 5 1 ^ ■ 



o 5 



0.5 



5 Or 



25 040 ^ 



T x o , 24 «1 a 4 

V = 2oiC -I — TO, 

o 5 



while to is a free parameter if 03 = 12 oi 04 / 05 or 



(03 05 — 2oi 04)^ 

TO= ^2V2 ■ 

300 a 4 a 5 

The shapes of the kink solutions (1.25), (1.27) are similar to that shown 
in Fig. 1.4, while the kink (1.28) may have also different profile shown by 
solid and dashed lines in Fig. 1.5. 
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Fig. 1.5 Kink-shaped waves with a ”hat” 



1 . 1.4 Periodic nonlinear waves 

Usually single bell-shaped solitary wave solutions are the particular cases of 
more general periodic solutions. Thus Korteweg and de Vries (1895) found 
the periodic solution of the KdV equation (1.1), 

u = 6^fc 2 ^1 — ft 2 + -p + n 2 cn 2 (kO , k)^ (1.29) 

where K and E are the complete elliptic integrals of the first and the 
second kind respectively , k is the modulus of the J acobian elliptic function 
Bateman and Erdelyi (1953-54); Byrd and Friedman (1954); Newille (1951). 
They called Eq.(1.29) the cnoidal wave solution since it is expressed through 
the Jacobi elliptic function cn . Cnoidal wave is not a linear superposition of 
the bell-shaped solitary waves. It tends to the single solitary wave solution 
(1.2) at k — > 1 as shown in the right column in Fig. 1.6 1 . Exact periodic 
and bell-slraped solitary wave solutions correspond in the same manner 
in case of the generalized 5th-order KdV equation (1.3) and the DMKdV 
equation (1.13). 

Although many equations, like Burgers’, BKdV and DMKdV equations 
have not exact bounded periodic solutions that transform into the kink- 
slraped ones, there exist exceptions. Thus the MKdV periodic solution 

1 Reprinted with permission from Elsevier Science 
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Ablowitz and Segur (1981) 



2d 



u= \ msn(m9,K), 



(1.30) 



transforms into the kink solution (1.27) at k — > 1. Another example is a 
dissipative nonlinear equation 



u t + 2 buu x + 3c u 2 u x + du xxx + f(u 2 ) xx + gu xx = 0, 



(1.31) 



This evolution equation represents an analog of the hyperbolic equation 
to be derived further in Sec. 5.3. Its bounded periodic solution Porubov 
(1996) is 



with 




cn(m 9, k) sn{rn. 9, k) dn{m. 9 , k) 
Ci + cn 2 {m9, n) 



b_ 

3c' 



Ci 



1 — 2k 2 + \/ k a — k 2 + 1 2 3 g 2 — V 

3k 2 ’ 4 V k 4 — k 2 + 1 ’ 



(1.32) 



and the following restrictions on the coefficients: 



/ 




b = 3 gy/—c. 



The periodic wave solution (1.32) has a functional form different from both 
the KdV cnoidal wave and the MKdV bounded periodic solution. When 
k = 1 we have Ci = 0, and the solution (1.31) tends to the kink-shaped 
solution (1.27) as it is shown in the left column in Fig. 1.6 in comparison 
with the transformation of the KdV cnoidal wave solution to the bell-slraped 
solitary wave. 



1.2 Formation of nonlinear waves of permanent shape from 
an arbitrary input 

All solutions presented in previous section require specific initial conditions. 
In practice more important is to know how an arbitrary finite amplitude 
input evolves. Analytical solutions of unsteady problems may be obtained 
if governing nonlinear equations are integrable Ablowitz and Segur (1981); 
Bhatuagar (1979); Calogero and Degasperis (1982); Dodd et. al (1982); 
Newell (1985), otherwise only numerical solutions are available. Often ini- 
tial input transforms into the stable quasistationary wave structures of 
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a 




Fig. 1.6 Comparison of the periodic solution (1.32) (left column) and the KdV cnoidal 
wave (1.29) (right column) for different values of the Jacobi modulus: (a)/-e 2 = 0.25, (b) 
k, 2 = 0.995, (c )k 2 = 0.99995, (d)/^ 2 = 1. After Porubov and Velarde (2002). 



permanent form which may be described by the analytical solutions. Also 
the analysis gives the conditions when the formation of them is possible. 
In this section we illustrate it using some instructive examples. 
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1.2.1 Bell-shaped solitary wave formation from an initial 
localized pulse 

The fifth-order KdV was extensively studied numerically. The oscil- 
latory travelling solitary- wave solutions were found in Boyd (1991); 
Nagashima and Kuwahara (1981). The evolution of the initial monotonic 
solitary wave into radiating or oscillatory solitary waves was simulated in 
Benilov et. al (1993); Karpman and Vanden-Broeck (1995). In a series of 
papers Salupere et. al (1997); Salupere et. al (2001) the solitary wave forma- 
tion from a periodic input was studied for an equation similar to Eq.(1.3). 
We shall study the evolution of a localized initial pulse. Previously, local- 
ized pulse evolution into an oscillatory solitary wave was considered in 
Nagashima and Kuwahara (1981) for the equation ut + uu x — y 2 u$ x = 0. 
Below we consider the formation of solitary waves in the systems governed 
by Eq.(1.3). Following Porubov et. al (2002) we use two methods for com- 
putations, finite-difference and pseudo-spectral, see Sec. 2.3.2. Below only 
those results are shown that were obtained using both numerical methods. 
We have tried various shapes of the initial localized pulses, rectangular, 
Gaussian distribution etc. 

Influence of the fifth-order dispersive term. First the 5th-order KdV 
equation was studied. Since the role of the fifth-order derivative term is of 
interest the coefficients b and d in Eq.(1.3) were fixed for all computations, 
b = 1, d = 0.5. We found that the rectangular initial pulse splits into a 
sequence of solitary waves when the coefficients of dispersive terms, d and 
/, are of opposite sign. For both coefficients positive the initial rectan- 
gular profile is dispersed without formation of any localized waves. The 
dependence upon the sign of the ratio d/ f is in agreement with the exact 
solitary- wave solution (1.4) and the analysis of the dispersion relation 
Karpman (1993). However, more smooth Gaussian initial profiles provide 
the appearance of solitary waves even for positive coefficients when / is 
rather small, e.g., / = 0.01. 

The next result we have obtained is the dependence of the number of 
solitary waves upon the value of / when d/ f < 0. 

Shown in Fig. 1.7 is the formation of the train of solitary waves from a 
Gaussian initial pulse in the KdV case, / = 0. 

Figs. 1.8-1.10 demonstrate the decrease of the solitary waves for / = 
— 1,-10,— 50 respectively. Both the amplitude and the velocity decrease 
with the increase of the absolute value of / in qualitative agreement with 
the exact solution (1.4). The ratio between the amplitude and the velocity 
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Fig. 1.7 Evolution of the initial Gaussian profile in the KdV case, / = 0. 



of each solitary wave in Fig. 1.7 is equal to 1.5 just as for the KdV soliton 
(1.2) . This ratio (and the amplitude) decreases with the decrease of /, from 
1.43 at / = — 1 to 1.33 at / = —50. A similar tendency is revealed by the 
phase analysis of single travelling wave solutions, cf. Fig. 2 in Kawahara 
(1972). The ratio for the 5th- order KdV exact solution (1.4) is 1.46, the 
amplitude and the velocity for given b and d are —105/(1352/), —9/(169/), 
respectively. Only at small / = —0.1 -4 — 0.15 are the numerical results for 
the leading solitary wave in quantitative agreement with the exact solution 
(1.4). More important is that the decrease of / affects the solitary- wave 
transformation from monotonic KdV solitons (1.2) to the oscillatory soli- 
tary waves. For convenience the last stages from Figs. 1.7-1.10 are collected 
in Fig. 1.11. In case / = —1, Fig. 1.11(B), the higher leading solitary wave 
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Fig. 1.8 Evolution of the initial Gaussian profile in the 5th-order KdV case, f = — 1. 



is oscillatory while other solitary waves remain monotonic, then the trans- 
formation occurs also for the second solitary wave, Fig. 1.11(C) . Alternate 
transformation of the solitary waves confirms the dependence of the kind 
of solitary wave upon the value of the product of / and the wave amplitude 
(hence, its velocity) found in Kawalrara (1972). 

Finally, we have found that simultaneous triggering of the signs of 6, d, 
f doesn’t affect the shapes of the solitary waves. They simply evolve to the 
opposite direction according to the analysis of the exact solution (1.4). 
Influence of the cubic nonlinearity. Now we add cubic nonlinearity, 
c / 0, while r = s = 0. The analytical solutions predict an action of 
the cubic nonlinear term depends upon the sign of c. Also the solution is 
sensitive to the ratio between nonlinear terms contributions, 6/c, and the 
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Fig. 1.9 Evolution of the initial Gaussian profile in the 5th-order KdV case, / = —10. 



value of the amplitude of an initial pulse. Indeed we have found that at 
b = 1, c = 1, d = 0.5, f — —1 the train of solitary waves arises only from 
positive initial pulse with amplitude equal to 0.5 while a negative one is 
dispersed. Dependence of the amplitude on the sign of b is very typical for 
the exact solutions of the equation with quadratic nonlinearity, see (1.2), 

(1.4) . However, at smaller b, b = 0.2 or b = —0.2, the formation of solitary 
waves no longer depends upon the sign of the initial pulse amplitude. A 
similar tendency is observed when the initial amplitude becomes close to 1 
or higher when predominant cubic nonlinearity excludes an influence of the 
quadratic one on the sign of the wave amplitude like in the exact solution 

(1.5) , (1.6). 

At the same time, the stage of forthcoming evolution of already gen- 
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Fig. 1.10 Evolution of the initial Gaussian profile in the 5th-order KdV case, / = —50. 



eratecl solitary waves is not so sensitive to the value of 6, while other an- 
alytical restrictions on the coefficients become more important. Besides 
the condition dj f < 0 following from the linear analysis Karpman (1993); 
Kawahara (1972), there is f/c < 0 resulted from the nonlinear exact so- 
lution (1.5), (1.6). Moreover, at small / one can anticipate an evidence of 
the condition d/c > 0 given by Eq.(1.8). 

When c > 0, / < 0 all above mentioned inequalities are satisfied. The 
number of solitary waves generated from the initial localized pulse increases 
with the increase of the value of c and fixed values ofd=0.5,/=— 1 and 
also b = 1 or b = 0.2. The velocity of the waves increases also, the width 
(proportional to 1 /k) decreases, while the amplitude remains practically one 
and the same. We also observed the alternate transformation of the solitary 
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A 




x 



B 




x 



C D 





Fig. 1.11 Transformation of the kind of solitary waves in the 5th-order KdV case. (A) 
/ = 0, (B) / = -1, (C) / = -10, (D) / = -50. 



waves from monotonic to oscillatory when c increases for both values of 6. 
Independence of the amplitude of c doesn’t follow from the exact solution 
(1.6) as well as from the asymptotic solution. The third formula in (1.6) 
predicts growth of positive values of k 2 only for b = 0.2 giving negative 
values for 6=1. However, let us express k through the amplitude A and 
substitute it into the expression for the velocity V. Then one can exhibit for 
both values of 6 the similarity of the variation of the velocity with respect 
to c with that obtained in numerics. 

As found in previous subsection, the decrease of the negative / values 
affects the decrease in the number of solitary waves. Assume 6 = 1, d = 0.5, 
we have tried simultaneous variations of c and / in order to sustain one 
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A 




x 



B C 





Fig. 1.12 Conservation of the number of solitary wave thanks to the simultaneous pres- 
ence of cubic nonlinearity and fifth-order derivative term. (A) c = 0, / = 0, (B) c = 15, 
/ = -10, (C) c = 75, / = -50. 

and the same number, see Fig. 1.12. We see that the wave amplitude 
keeps its value from Fig. 1.12(A) to Fig. 1.12(C) while the velocity growths. 
This confirms that the amplitude depends upon the ratio f /c but velocity 
is proportional to /. The kind of solitary wave alters from monotonic 
Fig. 1.12(A) to oscillatory Fig. 1.12(B, C). 

At c < 0, d > 0, / > 0 only one inequality, f/c < 0, is satisfied, and 
solitary wave formation is observed only for small / , otherwise the initial 
pulse is dispersed in this case. On the contrary small absolute values of c 
provides the solitary waves formation at c < 0, d > 0, / < 0. When both c 
and / are positive no solitary waves appear. 

Influence of nonlinearities s UjU^ and r uu lra . Suppose c = r = 0 
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r=-l . 50 r=-l . 57 





r=-l . 60 




X 



Fig. 1.13 Equalization of the first and the second solitary waves and subsequent ex- 
ceeding of the second wave due to the alteration of the negative values of the coefficient 
r. 



and vary s at fixed 6, d and /, that we choose b = 1, d = 0.5, / = — 1 . It 
is found that the amount of solitary waves and its transition from mono- 
tonic to oscillatory don’t depend upon the value of s. The wave amplitude 
decreases with the increase of s while the velocity keeps its value. Wave 
behavior is not sensitive to the sign of s. The condition for the solitary 
wave formation d/ f < 0 remains valid. 

The fact the velocity doesn’t depend upon s is in agreement with the 
exact solution. Certainly, solitary waves exist outside the restriction from 
(1.10). We also used numerical values of the amplitude to define k and 
then V using (1.10). A comparison of the velocities with those obtained 
numerically demonstrates the more agreement the less is the value of b. 
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Asymptotic solution (1.17), (1.19) also predicts the decrease of amplitude at 
permanent velocity. Indeed, we get that u max = tto(0)+5ui(0) = 6dk 2 /b(l+ 
k 2 [f/d — s/ (26)]) . At coefficient values we used the value of '« max decreases 
with the increase of s (s is not large in the asymptotic solution), while the 
exact solution predicts the same behavior only for positive values of s. 

When s = c = 0 the behavior of the solution differs from the previous 
one. Having the same values for b , d , / we obtain that increase of positive 
values of r yields a decrease in the velocity and an increase in the amplitude 
of the solitary waves. The number of solitary waves decreases. However, 
at negative values of r we found that at the initial stage of the splitting 
of the Gaussian profile the amplitude of the second solitary wave becomes 
equal to that of the first one at r = —1.57, see Fig. 1.13. At lesser r second 
solitary wave becomes higher, and two solitary waves form a two-humps 
localized structure shown in Fig. 1.14. 

It is no longer quasistationary since amplitudes of the humps vary in 
time. It looks like an interaction of two solitary waves when the second 
higher solitary wave surpasses the first one, then it becomes lower, and 
the process repeats. Decreasing r we achieve formation of a tlrree-humps 
localized structure shown in Fig. 1.15. Its evolution is similar to those 
presented in Fig. 1.14. Finally, only multi-humps localized structure arises 
from an initial pulse as shown in Fig. 1.16. The localized multi-humps 
structures in Figs. 1.14-1.16 keep their width, while their shapes vary in 
time. 

Certainly, unsteady multi-humps localized structures are not governed 
by the ODE reduction of Eq.(1.3) and, hence cannot be explained either by 
the phase portraits analysis or by the exact travelling wave solution (1.5), 
(1.9). Moreover, at negative values of / the exact solution doesn’t predict 
propagation to the right of the solitary wave with positive amplitude. 
Absence of linear dispersive terms. We have found an exact solitary 
wave solution that may be supported by higher -order nonlinear terms even 
without linear dispersive terms, at d = 0 or / = 0. Numerical simulations 
show that there are no solitary waves at both zero d and /. Some solutions 
from previous subsections keep their features at d = 0, in particular, this 
relates to the case r/ 0. At the same time cubic nonlinearity at d = r = 
s = 0 supports two- humps localized waves for c > 0. At negative c the wave 
picture is similar to those at d ^ 0. No stable solitary waves propagate in 
absence of only the fifth-derivative term, d ^ 0, / = 0 with the exception of 
the Gardner equation case where Slyunyaev and Pelinovsky (1999) found 
generation of the limiting amplitude solitons. 
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Fig. 1.14 Two-humps solitary wave formation at r = —1.6. 

To sum up, both higher order nonlinear and dispersive terms affect the 
formation of localized nonlinear waves their shape and their parameters. 
Thus, the number of solitary waves and the transition from monotonic to 
oscillatory wave are under responsibility of both 5tlr- order linear dispersive 
term, cubic nonlinearity and higher- order quadratic nonlinearity ru u xxx . 
More important is the formation of an unsteady but localized multi-humps 
wave structure thanks to r u u xxx and cubic nonlinearity at d = 0. The sign 
of the coefficient b of the KdV quadratic nonlinear term is important for 
choosing the sign of the input amplitude. At the same time the nonlinearity 
su x u xx doesn’t affect the formation and behavior of solitary waves. 

Certainly, the shapes of the resulting solitary waves are not obviously 
governed by the exact and asymptotic travelling wave solutions. Some 
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Fig. 1.15 Three- humps solitary wave formation at r = —3.2. 



other features of numerical solutions, like the dependence of the number of 
solitary waves upon the values of the equation coefficients or a transition 
from monotonic wave to an oscillatory one, are not predicted by analytical 
solutions. However, the combinations of equation coefficients required for 
the existence of solitary wave are realized in numerics. Also numerical 
wave amplitude and velocity relate like in the analysis. Evidence of all 
these predictions even qualitatively is very important for a justification of 
the numerical results. 

Formation of the bell-slraped solitary waves in presence of dissipation 
or an energy influx will be considered further in the book. 
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Fig. 1.16 No solitary waves other than multi- humps one at r = —12. 



1 . 2.2 Kink-shaped and periodic waves formation 

The formation of the kink-shaped waves was studied considering the 
evolution of the Taylor shock from discontinuous (step) initial condi- 
tions under the governance of the Burgers equation Saclrdev (1987); 
Wlritham (1974). It was found the appearance in due time the steady state 
kink solution (1.22). A quasihyperbolic analog of the Burgers equation was 
studied in Alexeyev (1999) where it was found that kink may be formed 
from suitable initial conditions. The Korteweg-de Vries-Burgers equation 
(1.24) also extensively studied but mainly in the direction of generation of 
the triangle profiles and oscillating wave packets Berezin (1987), see also 
references therein. Of special interest is the formation of the kinks with a 
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”hat” shown in Fig. 1.5. One possibility will be considered in Sec. 5.3. 

Usually periodic waves are generated in finite domains from a harmonic 
input. Thus the KdV cnoidal waves (1.29) are realized numerically and in 
experiments in a paper by Bridges (1986), also Kawalrara (1983) obtained 
numerically periodic wave structures in a system governed by Eq.(1.13) 
with 05= 0, while at nonzero coefficient similar results were found in 
Rednikov et. al (1995). Note that harmonic input in the finite domains is 
used also for the study of the bell-shaped solitary waves interactions where 
no periodic wave structure of permanent shape arises Salupere et. al (1994); 
Salupere et. al (1997); Salupere et. al (2001). 



1.3 Amplification, attenuation and selection of nonlinear 
waves 

As already noted the bell-shaped solitary wave is sustained by a balance 
between nonlinearity and dispersion. What happens with the wave when 
dissipation/accumulation destroys this balance? It was shown in Sec. 1.1 
that the bell-shaped wave may exist even in presence of dissipation but un- 
der strong restrictions on the equation coefficients. Assume the influence 
of dissipation/ accumulation is weak and is characterized by a small pa- 
rameter e << 1. It turns out that an asymptotic solution may be found in 
this case whose leading order part is defined as a solitary wave with slowly 
varying parameters. Depending on the problem either slow time, T = et, 
or slow coordinate, X = ex, may be used. 

In the former case the solitary wave solution is 

u{0, T ) = A(T) cosh ~ 2 {k{T)9), (1.33) 

where 9 X = 1, 9 t = —V(T). In the latter case we have 

u{9, X) = A(X) cosh ~ 2 (k(X)9), (1.34) 



where 9 X = P(X), 9 t = —1. Next order solutions give us the functional 
form of the dependence of the wave parameters upon the slow variable. 
Derivation of the asymptotic solution will be described in Sec. 2.2. Now only 
general features of the wave behaviour are considered. When A(T) increases 
in time k(T) usually increases also, hence the width of the wave inverse 
proportional to k(T), decreases. This is an amplification of the solitary 
wave. On the contrary, we have an attenuation of the solitary wave when its 
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Fig. 1.17 Temporal evolution of an initial solitary wave resulting in a selection: a) from 
below, b) from above. 



amplitude decreases while its width increases. Sometimes it happens that 
the increase/decrease of A(T) takes place not up to infinity/zero but to 
the finite value A*. Usually this value is defined by the governing equation 
coefficients, hence, by the physical parameters of the problem. To put 
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Fig. 1.18 Solitary wave (1.33) (dashed line) vs solitary wave (1-34) ( solid line). 

this another way, the parameters of the resulted steady wave are selected. 
Selection provided by an amplification of an initial wave, will be called 
selection from below, see Fig. 1.17(a), while selection from above happens 
as a result of an attenuation of an initial wave, see Fig. 1.17(b). 

Shown in Fig. 1.18 is the profile of the wave (1-34) in comparison with 
the symmetric solitary wave solution (1.33) at t = 0. One can see that the 
wave (1-34) is asymmetric with respect to its core (or maximum). However, 
only initial stages of the temporal evolution of (1-34) differs from that of 
Eq.(1.33). As follows from Fig. 1.19, the final stage of the selection both 
from below and above, is the symmetric bell-shaped solitary wave like shown 
in the last stages in Fig. 1.17. 

The amplification/attenuation of the kink may be described asymp- 
totically and numerically Sachdev (1987), it will be shown in Sec. 5.3. 
Cnoidal wave evolution may be accounted for an asymptotic solution 
similar to that of the bell-slrapecl solitary waves Rednikov et. al (1995); 
Svendsen and Bulrr-Hansen (1978). 
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Chapter 2 

Mathematical tools for the governing 
equations analysis 



As a rule governing equations for nonlinear strain waves are nonitegrable 
by the inverse scattering transform method, and only particular analytical 
solutions may be obtained. Hence the study of real physical processes 
requires a combined analytico-numerical approach. The aim of this chapter 
is to describe methods to be used in this book. The choice of the analytical 
and numerical procedures is based on an experience of the author and does 
not claim a completeness. 



2.1 Exact solutions 

2.1.1 Direct methods and elliptic functions 

Most of the mathematical work in the realm of nonlinear phenomena 
refers to integrable equations and their exact solutions, particularly, pe- 
riodic. Among the recently developed general methods the algebroge- 
ometrical approach may be used in an efficient way to find such solu- 
tions. Not only the numerical realization and graphical representation of 
the solution is provided by this method but also multi phase quasiperi- 
odic solutions as well as purely periodic ones may be represented using 
the algebrogeometrical approach as illustrated in Belokolos et. al (1994). 
When we are interested in a self-similar solution of a partial differential 
equation one can use well developed theory of the solutions of ordinary 
differential equations, see, e.g., Ince (1964). Exact solutions of nonlin- 
ear nonintegrable partial differential equations are obtained usually us- 
ing various direct methods. The significant point in direct methods is to 
build in advance the appropriate functional form of the solution (ansatz) 
of the equation studied. For example, the usage of ansatz in the form 
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32 Amplification of Nonlinear Strain Waves in Solids 

of a hyperbolic tangent (tanh) power series resulted in finding of new 
exact travelling wave solutions (see, e.g., Korpel and Banerjee (1984); 
Parkes and Duffy (1996) and references therein). The choice of tanh 
is caused by the fact that any derivatives of tanh may be expressed as 
a polynomial with respect to the tanh itself. Then the equation stud- 
ied becomes polynomial of the tanh after substituting the ansatz, and 
solution parameters are obtained from the algebraic equations appear- 
ing after equating to zero coefficients at each order of tanh . One would 
like to apply the same procedure to find more general periodic solutions. 
First of all another appropriate ’’basic” function is required instead of the 
tanh . For this purpose various elliptic functions were proposed recently, 
the most popular were theta functions (see, e.g., Chow (1995), Nakamura 
(1979)), Jacobian elliptic functions (see, e.g.,Kostov and Uzunov (1992); 
Parker and Tsoy (1999)) and the Weierstrass elliptic function Kasclreev 
(1990); Porubov (1993); Porubov (1996); Porubov and Parker (2002); 
Samsonov (1995); Samsonov (2001). In principle, periodical solutions could 
be obtained in terms of any of these functions. The efficient choice is caused 
by the simplest procedure of the ansatz construction and the least compli- 
cated algebra for determining solution parameters. It is well known that 
theta functions may be included in the Hirota bilinear method in order to 
get N-periodical solutions Nakamura (1979). However most of dissipative 
equations cannot be transformed to the bilinear form. At the same time 
single travelling wave solution derivation looks very complicated even for 
non-dissipative equations Chow (1995), Nakamura (1979). Moreover we 
have to deal with four theta functions that results in additional difficul- 
ties for the ansatz construction. Explicit periodic travelling wave solutions 
may be found for many nonintegrable equations and systems by using an 
ansatz in terms of p, with appropriate forms for the ansatz suggested by 
information about the poles of the solution. When compared to the use 
of theta functions or Jacobi elliptic functions, a prime advantage of using 
the function p is that the algebra is drastically simplified. The ansatz for 
the solution involves only one Weierstrass function p(C, < 72 , < 73 ), instead of 
four theta functions or three Jacobi elliptic functions cn(£, fc), sn(f, k) 
and dn(C, k). Another advantage is that two apparently distinct solutions 
are readily recognized as equivalent. In order to see it let us first give 
some properties of the Weierstrass elliptic function p to be used below. 
According to its definition Whittaker and Watson (1927), the Weierstrass 
function is analytical in the complex plane other than in the points where 
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Mathematical tools for the governing equations analysis 33 

it has double poles. The governing equation for the function p is: 

{p , (C)} 2 = 4p 3 - 52 p-ff3 (2.1) 

where <72 and 93 are constant parameters. Remarkable features of the func- 
tion p are that all of its derivatives can be written by means of itself, 
and that any elliptic function f may be expressed using p it and its first 
derivative as Whittaker and Watson (1927) 

/ = A(p) + B{p)p', (2.2) 

where A and B are rational functions with respect to p. Depending on 
the ratio between <72 and <73 the Weierstrass function may be bounded or 
unbounded inside the domain under study. The bounded periodic solutions 
are more conveniently expressed by writing them in terms of the Jacobi 
elliptic functions cn, sn and dn which are bounded on the real axis. For 
this purpose the relation between the Weierstrass function and the Jacobian 
functions is used as a special case of (2.2). Indeed, the familiar link is 

obtained in Whittaker and Watson (1927) but using the singular function 

_2 
sn , 

P«> 92, 93) = e 3 + (ei - e 3 )s7W 2 (v / e 1 - e 3 C, k). (2.3) 

However, following the method introduced in Whittaker and Watson (1927) 
one can check that the following formula is valid: 



p(C. 92, 93) = e 2 - (e 2 - e 3 )cn 2 (v'e 1 - e 3 (, k ), (2.4) 

connecting the Weierstrass function with the Jacobi function cn, regular 
along the real axes. Here k = \/{e 2 — e 3 )/(ei — e 3 ) is the modulus of the 
Jacobian elliptic function, while t = e m (m = 1, 2, 3 , e 3 < e 2 < e± ) are 
the real roots of the cubic equation 



4r 3 - g 2 r - g 3 = 0. 



(2.5) 



Expressing these results in terms of an appropriate choice of parameters, 
the wave number k = y/ei — e 3 and the Jacobian elliptic modulus k, we 
have 



e 3 = 



1 + fc 2 



2 k 2 - 1 



2 — k 2 



e2 = 



ei = 
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92 = !|k 4 (1 - k 2 + fc 4 ), g 3 = ^« 6 (fc 2 + 1 )( 2 - k' 2 ){l - 2k 2 ). (2.6) 

The localized both the bells-shaped and the kink-shaped solitary wave so- 
lutions appear in the limit k — » 1 of the Jacobi elliptic functions. 

Now consider an instructive example Porubov and Parker (2002). In 
Parker and Tsoy (1999), solutions were sought in terms of powers and 
products of Jacobi functions and thereby two solutions were obtained 

zi= e 2 — (e 2 - e 3 ) cn 2 ( v / e i - e 3 C, k), (2.7) 



r 2 m 2 2 r 2 m 

Z 2 = — - — sn, yrQ,m) ± — — cn{rQ, m)an{rQ, m) — 



r 2 (l + rn 2 ) 



12 



( 2 . 8 ) 



which appear different. Obviously, the solution (2.7) is a representation of 
the Weierstrass function (2.3). However, one can check by direct substitu- 
tion that the solution (2.8) also satisfies equation (2.1) when the parameters 
r and to are defined as solutions to 

^(1 + 14to 2 + to 4 ) = 92 , ^g(i - 33 to 2 - 33to 4 + to 6 ) = g 3 . (2.9) 

Therefore z 2 is also a solution satisfying the same governing equation 
(2.1) defining the Weierstrass function. The two expressions (2.7) and (2.8) 
are essentially equivalent solutions, provided that (f, to) are appropriately 
related to (r, k) . Accordingly, working with p reveals links between seem- 
ingly distinct forms of solution. 

It is to be noted that the first Weierstrass function derivative p' can- 
not be expressed as a polynomial of the Weierstrass function itself Whit- 
taker and Watson (1927), and we have to equate zero separately coeffi- 
cients at each order of p and at products of p' and corresponding orders 
of p Kascheev (1990); Porubov (1993); Porubov (1996); Samsonov (1995); 
Samsonov (2001). Therefore we really deal with two functions, and it is 
unlikely to get the solution using the ansatz proposed in the form of power 
series with respect both of the p and p' as it was done for the tanh. An- 
other idea may be used, based on the singular point analysis of the possible 
solution and the well known fact that the Weierstrass function p has the 
second order pole Whittaker and Watson (1927). In order to check the poles 
of a solution we shall use the WTK method Weiss et. al (1983) looking for 
the solution in the form of Laurent- type series about the singular manifold. 
It will be explained in details in Sec. 2.1.2. 
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Mathematical tools for the governing equations analysis 35 

The exact solutions obtained in this manner, usually belong to the class 
of travelling wave solutions which require special initial conditions. In case 
of the solitary wave solution the initial condition should be have the shape 
of the solitary wave itself. Moreover, travelling wave solutions for the dis- 
sipative equations usually have not free parameters, and additional rela- 
tionships on the equation coefficients are required for the existence of the 
solutions Porubov (1993); Porubov (1996); Porubov and Velarde (1999); 
Porubov and Parker (1999); Porubov and Parker (2002); Samsonov (1995); 
Samsonov (2001). 



2 . 1.2 Painleve analysis 

Recently it was developed the theory of nonlinear ordinary differential 
equations whose solutions have not movable singularities, other than poles. 
Then the theory has been extended to partial differential equations. Usually 
these equations are called ’’equations with the Painleve property”. The 
achievements of the theory may be found in Cariello and Tabor (1989); 
Conte (1989); Conte et. al (1993); Levi and Winternitz (1992); Newell et. al 
(1987); Weiss et. al (1983). Here we concentrate on the one aspect of the 
tlreory-tlre singular manifold method or WTC method for partial differen- 
tial equations Newell et. al (1987); Weiss et. al (1983). Let tp(x,t) = 0 is 
the ’’singular” or ’’pole” manifold on which a solution u(x,t) is singular. 
The main idea of the WTC is to demonstrate that the expansion 

u(x,t) = — (2N0) 

^ j=o 

is single valued. This requires (i) a is an integer, (ii) is analytic in x 
and t and (iii) the equations for the coefficients qj have self-consistent 
solutions. When all these conditions are satisfied the equation under study 
has the Laurent property. Also it is necessary to assume that neither ip x 
nor ip t vanish on ip(x,t) = 0. Let us illustrate how the methods works 
on an example of the KdV equation (1.1). Following Newell et. al (1987); 
Weiss et. al (1983) we assume b = 3, d = 1. Substituting the ansatz (2.10) 
into Eq.(l.l) one can find a = 2 and q 0 — — 2(^. Recursion relations for 
the q.j are 

(j + !)(i - 4 )(j - 6 )qj = F(ip t ,<p x , ..., q k ; k < j). (2.11) 

The values of j , j = — 1, 4 and 6 are called resonances. At each such 




February 11, 2004 20:10 



WSPC/Book Trim Size for 9in x 6in 



ws-book9x6 



36 Amplification of Nonlinear Strain Waves in Solids 

resonance the right-hand side of Eq.(2.11) vanishes thus ensuring the inde- 
terminancy of the corresponding qj. Moreover, the expansion (2.10) may 
be truncated at O(<p 0 ). As a result we obtain using Eq.(2.11) an auto- 
Backlund transformation for the solution of Eq.(l.l), 

. „ d , 

u{x,t) = 2 —\ogp + q 2 , 

where q 2 satisfies the KdV equation. Also the Lax pair for the KdV equation 
follows from the solution of Eq.(2.11) Newell et. al (1987); Weiss et. al 
(1983). 

However, the Lax pair cannot be obtained for nonitegrable equations as 
opposed to a truncated expansion that carries an information about possible 
pole orders of a solution. Using this information the anzats for the solution 
may be proposed, in particular, in terms of the Weierstrass elliptic function 
p and p'. Substituting the proposed form of the solution into the equation 
under study and equating to zero coefficients at each order of p and at prod- 
ucts of p' and corresponding orders of p one can get the algebraic equations 
on the solution parameters, the phase velocity and the Weierstrass function 
parameters g 2 ,gs ■ Certainly this procedure is of phenomenological kind 
but it allows to obtain the solutions of nonitegrable nonlinear equations in 
an explicit form. Some examples are presented below. 

2.1.3 Single travelling wave solutions 

First we consider exact solutions of DMKdV Eq.(1.13) obtained in Porubov 
(1993). It was found there the following auto-Backlund transformation for 
the solution u(x, t) 

12 «4 . 12 . 2040 : 4 .., , ~ . 

u= — — (log (p)xx H —((*3 ~ )(log<p) x + u, ( 2 . 12 ) 

05 5 o 5 05 

where u (x,t) satisfies Eq.(1.13). Concerning only tavelling wave solutions, 
one can reduce Eq.(1.13) to the third-order ODE of the form: 

a. 4 u" + 03 u"+ a 2 u'+ 05 uu + 01 u 2 — V u + P = 0, (2.13) 

where P = const , prime denotes differentiation with respect to 9,9 = x — 
Vt. Based on Eq.(2.12) possible solution may contain simple and second- 
order poles that may be modelled in terms of p as Porubov (1993): 

Bp' 



u = Ap + 



P + C 



+ D. 



(2.14) 
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Substituting Eq.(2.14) into Eq.(2.13) one can derive a system of alge- 
braic equations in A , B , C, D , phase velocity V and Weierstrass function 
parameters g 2 , g 3: 

(32C - <73 - 4 C 3 )B = 0, (12C 2 - g 2 )B = 0, 

P = VD + 8ai B 2 C - a\D 2 - a 3 g 2 A/2 - 2 a 2 BC + 12 a 4 BC 2 - 
a 5 (2 ABC 2 + g 2 AB/2 - 2 BCD), 

2ai(2 B 2 + AD) - VA + 2 a 2 B + 2 a 5 ( BD - ABC) = 0, 
oz±A 2 T 6a 3 ^4 12 a 4 B T 6 as AB — 0, 



B{2a\{D - AC) - V) = 0, A( 12 a 4 + a 5 A) = 0, 



2cx 4 AB 2a 3 -B+ a 2 /I — a 5 ( AD -t- 2 as B~) — 0. 

The solutions of these equations are: 

(i) when g 2 , g 3 are free parameters and a 3 = 2a 4 a 4 / as; 

A = B = 0,D = -b, V = - 2ai ° 2 . 



as 



«5 



«5 



(ii) when either 



C = — 



1 



r(a 3 — 



2a 1 a 4 . 2 
~ ) ) 



300 a 4 

or a 3 = 12a 4 a 4 / as, C is a free parameter 



12 a 4 6 2ai a 4 

A = — , B = — ( a 3 ), 



a 5 



5 a 5 



a 5 



1 



a 2 2a 4 2a 4 a 4 

D — — — + yyy (a 3 ~ ) + 

as as as 25 a 4 as 

24a 1 a 4 2a 4 a 2 Aa\ 2a 4 a 4 
V = — C- — 1- yy2-(a: 3 — ) + 



2a 4 a 4 2 

(as ~ ) , 



a 5 



«5 

2a 



«5 



a 5 



as 



25 a 4 a5 



^(as-^ 1 ) 2 , 32 = 12C 2 , g 3 = 8C 3 . 



a 5 
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38 Amplification of Nonlinear Strain Waves in Solids 

Using (2.4) the solution (2.14) with parameters defined by (i) may de- 
scribe a particular bounded cnoidal wave, propagating with fixed phase 
velocity, of the form: 

u = k 2 K 2 cn 2 (k6 , k) — — -ffk 2 (2k 2 — 1) (2.15) 

«5 CC 5 «5 

where k is a free parameter. When the Jacobian elliptic functions modulus 
k — > 1 the solution (2.15) transforms to the solitary wave solution (1.12). 
The solution (2.14) with parameters defined by (ii) accounts for a bounded 
kink-shaped solution (1.28). When C is a free parameter, kink propagates 
with any phase velocity value. 

Besides bounded solutions (2.15) and (1.28), the solution (2.14) also 
describes unbounded ones in the form of localized and periodic discontinu- 
ities. Finally it is to be noted that the functional form (2.14) in terms of 
the Weierstrass function is not unique. One can see that there exist at least 
one more solution of the form: 

u= ^r(a 3 - 2a l ) 2 exp(2y)p(exp(y) + G, 0, g 3 ), 

25 0 : 40:5 05 

where 

y = exp ( 76 O, 7 = xH a 3 - ^ ° 4 ) 1 

5 0:4 a 5 

G and <73 are free parameters. This solution allows to describe only the 
bounded kink-shaped wave (1.28). At the same time it accounts for a new 
periodically discontinious solution. 

In previous Chapter the bounded periodic solution (1.32) of the equa- 
tion (1.31) was considered. When studying travelling wave solutions, i.e. 
solutions depending only on the phase variable 6 = x — Vt, this equation 
may be transformed into the O.D.E., which results in the following equation 
after integrating once with respect to 9: 

dr]"' + gv" ~ Vri' + b(r] 2 y + f(r] 2 )" + c(r] 3 y = N. (2.16) 

Further we follow the results obtained in Porubov (1996). The transforma- 
tion is obtained here using the WTK method of the form 

u = f ± ^ f 2 ~ 2cd (log ip)' + u. (2.17) 

c 
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Therefore possible solution should contain first order pole. The Weierstrass 
elliptic function p possesses second order pole, and we shall propose three 
solution forms. The first of them is 

u = + (2.18) 

In order to find the solution parameters the formula (2.18) is substituted 
into the Eq.(2.16). Then equating to zero coefficients at each order of p+C 
and p' one can derive the algebraic equations on A , B , C, phase velocity V 
and the Weierstrass function parameters g 2l 53: 

(p + C)- 4 :(g 2 C-g 3 -4C 3 ) 2 (6Af-3A 2 c-6d) = 0, 

(p + C)~ 3 : ( g 2 C -g 3 - 4C 3 )(12C 2 - g- 2 )(UAf - 9A 2 c - 12 d) = 0, 

(p + C)~ 2 : 4 (g 2 C - g 3 - 4C 3 )(2bB — V + 48 AC f + 3c(12A 2 C - 
B 2 ) — 48Cd) — 3(12C 2 — g 2 ) 2 (2f — 2Ac — cL) = 0, 

(p + C)~ l : 4 (g 2 C -g 3 - 4C 3 )(10/ - 4/ + 3Ac + 6 d) - (12 C 2 - 

g 2 ){2Bb —V + 24 AC} + 3 {B 2 - 12 A 2 C)c - 12 Ad) = 0, 
p° : 4(12C 2 - g 2 )qA 2 = N, 

(p + C) : 2Bb - V + 24 AC f + 6 {B 2 - 12^ 2 C)c - 12 Cd = 0, 

(p + C ) 2 : 2 Af + d+ 2 A 2 c = 0 , 
p'(p + C)~ 3 : ( g 2 C -g 3 - 4 C 3 ){g + 2 Bf - Ab - 3 ABc) = 0, 
p'(p + C )~ 2 : (12C 2 - g 2 ){2Ab -g- 2 Bf + 6 ABc) = 0, 
p’ : 2 Ab + g + 2 Bf + 6ABc = 0. 

Three solutions of algebraic equations are obtained. The first appears 
when 

12C 2 - <72 = 0, g 2 C -g 3 - 4 C 3 = 0. 

In this case two of three roots e, of Eq.(2.5) are equal to one another, and 
no periodical solution exists. Then the solution (2.18) will have the form 
of localized discontinuity under positive C values. When the parameter C 
is negative we get k = 1, and the bounded kink-slraped solitary wave 

solution follows from Eq.(2.18): 

u = 7tanh(m0) + u 3 . 

For the wave amplitude 7 two formulaes are valid 

7 = Ann , A 3 = (/ + (/- 2cd) 1/2 )/(6c), 



(2.19) 
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7 = A 2 m, A 2 = (f — (f — 2cd) 1/2 )/(6c). 
Then for u o we have 

2 bAj — g 

U ° = 2(3 cAj + f) ’ J = ’ ’ 
and phase velocities Vi, V2 are 

Vl = 2uqi6 + 60 (^ 0 ! + 4 Aim 2 ) — 4m 2 (2/Ai — d), 



V 2 = 2u 0 2b + 6 c ( wq 2 + 4A%m 2 ) — 4m 2 (2fA 2 — d ), 

while to 2 = —3 C = 3ei is a free parameter. 

The second solution corresponds to the situation when 

g 2 C - g :i - 4C 3 = 0, 12C 2 - g 2 ± 0. (2.20) 



In this case the solution may exist under additional conditions on the equa- 
tion coefficients / = g = 0, that results in absence of the dissipative terms 
in Eq.(2.16). Hence it becomes now the O.D.E. reduction of the Gardner 
equation. The bounded cnoidal wave solution arises when C = —e\ and 
has the form 



u 



± 




cn(m9 , K)sn(md, k) 
dn(m9 , k) 



b_ 

3c’ 



( 2 . 21 ) 



where in 2 = ei — e^. It governs the travelling cnoidal wave, propagating 
with the phase velocity V = —b 2 /(3c) — 6eid that transforms into the 
kink-slraped soluiton (2.19) when k = 1. 

Finally, the third solution arises when 



12C 2 - ff2 = 0, g 2 C - g :i - 4C 3 ± 0. 



In this case the bounded cnoidal wave solution (1.32) holds. When k = 1 
we have C* = 0, and it transforms into the kink-slraped solution (2.19), see 
Fig. 1.6. 

Now we shall consider the second possible solution’s form in terms of 
the Weierstrass elliptic function p: 

g = \J Ap + B. (2.22) 

Substitution Eq.(2.22) into the Eq.(2.16) allows us to conclude that solution 
may now exist only when / = g = 0. Then one can get the algebraic 
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equations for the solution parameters equating to zero coefficients at the 
terms p' p k , k = 0 + 3, and ( Ap + B ) 5 / 2 : 

p' p 3 : cA + 2d = 0 
p' p 2 : 9cAB -VA + +15 dB = 0 
p' p : 9cAB - 2VA + +12 dB = 0 
p' : 12 cB 3 - 4VB 2 + dA(Bg 2 - Ag 3 ) = 0 
(. Ap + B) 5 ' 2 : N = 0 



One can obtain the following solution of these equations: 



A = 




2 dei 

, t 

c 



i + 3,y = 




where e, are the real roots of Eq. (2.5). We again deal with the O.D.E. 
reduction of the Gardner equation, however now one can obtain its another 
bounded solutions. Thus for d/fc> 0 we find cnoidal wave solutions of 
the form: 

(I) for B = 2de 2 /c 



u 




to k cn{m9 , k), 



(2.23) 



where to 2 = eq — e 3 and for the phase velocity we have + = dm 2 (2k 2 — 1). 
(II) for B = 2de\jc we obtain 



u 




to k dn(m6, k), 



(2.24) 



Solution (2.24) represents cnoidal wave propagating with the phase velocity 
V = dm 2 ( 2 — k 2 ). It is to be noted that at k = 1 both the solutions 
(2.23) and (2.24) transform to the bell-shaped solitary wave 



u = \ — to ch 1 (to0, k), 



propagating with the velocity V = dm 2 . 

When d/c > 0 the bounded solution appears from (2.22) if B = 2de 3 /c 
and has the form 




to k sn(mQ,n), 



u 



(2.25) 
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For the phase velocity we find V = -Jm 2 (l + «; 2 ). In contrast to solutions 
(2.23), (2.24) the solution (2.25) transforms to the kink-shaped solution 
(2.19) when k = 1. 

Finally one can construct the third possible solution in terms of the 
Weierstrass elliptic function p depending now on new independent variable 
y = exp (k6). Changing variable in auto-Backlund transformation (2.17) 
one can see that the following solution may be proposed: 

u = A^-+B, (2.26) 

P+C 

Substituting it into Eq.(2.16) one can find that solution may exist under 
/ = 0. Therefore we now deal with the O.D.E. reduction of the Gardner 
equation with linear dissipation. The solution of the form (2.26) for the 
last equation was already obtained previously in Samsonov (1995). It was 
found that this solution may correspond either to discontinuous periodical 
solution or to the bounded kink-shaped solution of the form similar to the 
solution (2.19). 



2 . 1.4 Exact solutions of more complicated form 

The procedure of the obtaining exact solutions is based on the reduction 
to the ODE, hence only travelling wave solutions may be found. However, 
sometimes this procedure is applied to obtain only the part of a solution, 
and more complicated profiles are decsribed. Consider some examples. 
Recently the coupled nonlinear Schrodinger equations (CNLS) 

iW t + sW xx + (rjWW* + oUU*)W = 0, 
ipU t + rU xx + (?? -'UU* + aWW*)U = 0, (2.27) 

with r,s = ±1; p,<J,r] real, have attracted considerable interest because 
of their role in governing various physical wave-guiding systems. These 
equations are not generally integrable by the inverse scattering transform 
method, except in cases of high symmetry (e.g. p = l,r] = a = ±1). The 
machinery of the Weierstrass function is applied in the search for travelling 
wave solutions of equations (2.27) in the form 

IT = w(C)e iC , U = u{Qe i4> , (2.28) 

where 9 = x — ct, £ = £(0, t), (f> = (j>(0, t), with u, w, ( and <fi real. Substitu- 
tion (2.28) into equations (2.27) allows to separate the real and imaginary 
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parts. The solutions for u and w are sought in the form 

w = sjFp 2 + Ap+B, u = sjGp 2 + Dp + E, (2.29) 
while for £ and <f> we have 



C = \sc8 — sC\ J ^ + Co (t), 

4> = ±rcp 0 - rC 2 J ^ + 4> o(t) , (2.30) 

where C'i (i = 1,2) are free parameters. The most notable difference from 
previous solutions is the dependence of the “frequency” and “phase” on w 
and u, when Ci ^ 0. This introduces significant features into the profiles 
of the real and imaginary parts of W and of U . Thus, we find for W of the 
form: 



W = w exp i Y, 

where w is defined by Porubov and Parker (1999) 



w = \ A 



2 k 2 - 1 



K 2 — k 2 n 2 cn 2 (n9 , k) + B/A 



while for the phase function Y we have 

csO sCAl[p, n, k] 



2 






(2.31) 



(2.32) 



(2.33) 



II[y>, n, k\ is the elliptic integral of the third kind, ip = arcsin (sn(«0)), n = 
— K 2 k 2 /(B/A + e 3 ) and C'i is defined by 



A = 



2{<r rrj — s) 2 _ A 2 



77(1 - cr 2 ) 



B 



> C i ~ T 1 4 T - 92 -r + 53 



A 



B 



A 



(2.34) 



cr(E A — B D) — 3s B 2(sa — rr]) 

7= A ’ °= 1-^ ’ 

where B, E, g 2 and 7/3 are free parameters. The shape of the wave (2.31) 
depends strongly on the values of the parameters k, k and H\. An example 
of the behavior of Re W is shown in Fig. 2. 1(a) for the case Hi G [e 2 , 63], c ~ 
k and at t = 0. We get an interesting wave, consisting of the carrier wave 
slightly modulated and with superposed periodic disturbances. Therefore, 
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Re W 





Fig. 2.1 Exact solutions of the CNLS. a) profile with moving disturbances; b) almost 
envelope wave solution. 

the wave shape is not determined solely by the amplitude wave shape. 
When c considerably exceeds the wave number k we get the profile closer 
to the usual envelope wave solution, see Fig. 2.1(b). 

Another interesting profiles correspond to the exact solutions of the 
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complex Ginzburg-Landau equation, CGLE, 



nit + pu xx + q \u\ 2 u = u (2.35) 

where the constant coefficients arep = p r +i pi, q = q r +iqi, with pj , qj ^ 0, 
( u,p,q ) £ C, 7 £ R. The subscripts t and x denote temporal and spatial 
derivatives, respectively. This equation appears in the description of a large 
variety of physical phenomena, e.g., in nonlinear optics, non-equilibrium 
pattern formation, lasers, superconductivity etc. Like in case of CNLS 
equations, first of all, we decompose the solution u(x,t) in its amplitude, 
y , and phase, 6 , both real, 

u=y(C)e ie , (2.36) 

where £ = x — ct, 9 = Substituting (2.36) into (2.35) and equating 

to zero the real and imaginary parts one obtains two coupled equations for 
the functions y and z = 0q. Then periodic and pulse solutions of the CGLE 
may be found Porubov and Velarde (1999). The first derivative of the real 
part of the periodic solution with respect to x, v = ( Reu ) x , is 

v = — sin (o + arcsin(^-^ 1 ) N ) , (2.37) 

y \ Z ) 

with 

y = J— {fc?dn 2 (fciC, «) - <5i}, z = Vy 2 yl + c 2 - ( 2 - 38 ) 

V Pr 

When C — 0, the zeroes of the first derivative (2.37) are defined by 
the zeroes of the function y x and correspond to the zeroes of the Jacobi 
functions cn and sn. Their positions do not change in time, and for u we 
have harmonic temporal oscillations of the spatially periodic state defined 
by the amplitude parts of the solutions (2.36). However, the situation 
changes dramatically when C/ 0. In this case Z never vanishes, and the 
zeroes of the first derivative are defined by the zeroes of the sin function 
only, whose position vary in time. A typical situation is shown in Fig. 2.2. 
Fig. 2.2(a) shows a structure with four spatially repeated parts. During half 
of the time period the shapes of these parts vary, and we get in Fig. 2.2(d) 
a profile which is practically the mirror image to Fig. 2.2(a). Qualitatively 
this evolution does not depend upon the value of the modulus k of Jacobi 
functions. 
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Re u 




Re u b 




Re u 




Re u 




Fig. 2.2 Evolution of the periodic solution (2.36) Reu vs x for times t = mp^7r/(3p r 7), 
0 < m < 3 with k = 0.9, £i = 0.5. a)m = 0, b)m = 1, c )m = 2, d)m = 3. 
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Consider now the pulse solution Porubov and Velarde (1999), 



u = 



QAh 



k cosh 1 (/c^) expz0, 



(2.39) 



with 



k 2 = -; 



7 



h=p 2 r+ Pi, h — Pr Qi Pi q r - 



2 (Pi + 4 A p r - 4 A 2 Pi) 

The first derivative for the real part of (2.39) is 



( Reu) x — ^ 6Ali(l + 4 A ) tanh(fcC)sin ^0 — arctan(^)^ . 

(2.40) 

Thus, from (2.40) it follows that additional zeroes of the first derivative 
may appear if 



k > 



h ( arctan(l/(2A))\ 

6 Ah 6XP V 2 A ) ' 



(2.41) 



The evolution of the real part of the solution (2.39) is illustrated in Fig. 2.3. 
Again we see that two initial maxima in Fig. 2.3(a) disappear, Fig. 2.3(e), 
then an initial minimum at £ = 0 is changed into a maximum, while two 
minima arise, Fig. 2.3(f-h). Therefore, our solution is breather- like. If 
(2.41) is not satisfied, there is a pulse solution whose spatial behavior is 
determined by the function cosh _1 (fc£) only with one extremum at £ = 0. 



2.2 Asymptotic solutions 

Particular exact solutions are insufficient for understanding physical pro- 
cesses. One can see that many equations consist in generalisations of the 
integrable equations like the KdV equation. Sometimes additional terms 
may be considered as small perturbations. Often straight asymptotic ex- 
pansions are incorrect Cole (1968); Jeffrey and Kawalrara (1982); Nayfeh 
(1973), and the matching asymptotic procedure Ablowitz and Segur (1981); 
Kodama and Ablowitz (1981) may be applied to find a solution in this 
case. The basic idea of the perturbation approach is to look for a solu- 
tion of a perturbed nonlinear equation in terms of certain natural fast and 
slow variables. Often there is a need in only one fast variable, such as 
phase variable 0 in the unperturbed problem. Depending upon the prob- 
lem either slow time, T — et, or slow coordinate, X = ex, e « 1 are 
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Re u a Re u b 




Fig. 2.3 Evolution of the pulse solution (2.39) Reu vs x for times t = m7r/(3 , 0 < m < 8. 
a)m = 0, b)m = 1, c)m = 2, d)m = 4, e)m = 5, f)m = 6, g)m = 7, h)m = 8. 

introduced. Fast variable is generalized in a perturbed problem assuming 
either Q x = 1, 9 t = —V(T) or 6 X = P{X), 9 t = — 1. The functions V(T) 
or P(X) are defined to remove secular terms. Some of the appropriate sec- 
ularity conditions are formed from Green’s identity as follows. Assume a 
solution u = u(0, T) is of the form 
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u(9,T)= uo(0,T) + eui(Q,T) + ... (2.42) 

The leading order problem is nonlinear equation for w 0 whose solution 
is known. Then in the next order an inhomogeneous linear equation for iq 
holds 



L(ui) = F(uo). 

Here L(u\) is a linearized leading-order equation operator. Denoting by 
Vi(i = 1, ... ,M ) the M solutions of the homogeneous adjoint problem, 

L A { Vi ) = 0, i= 

where L A is the adjoint operator to L, we obtain that 
L(ui)vi - L A (vi)ui = Fvi 



is always a divergence. Using the boundary conditions it may be integrated 
to give the secularity conditions. In particular, when u, — > 0 as \9\ — > oo 
while Ui is bounded the secularity condition is: 



iFd9 = 0 . 



(2.43) 



Then either V(T) or P(X) are obtained from Eq.(2.43). Method may 
be used for finding perturbed cnoidal wave solutions. In this case in- 
tegration in (2.43) is carrying over the period Rednikov et. al (1995); 
Svendsen and Bulrr-Hansen (1978). It is to be noted that a simple qua- 
sistationary expansion is not uniformly valid on |x| — > oo. Complete 
solutions is obtained using matching quasistationary solution to a non- 
stationary one for large \9\, e.g. \9\ ~ 1/e Ablowitz and Segur (1981); 
Kodama and Ablowitz (1981). 

Let us apply this method for a perturbed solitary wave solution of the 
DMKdV equation (1.13). Following Kliaklrandler et. al (2000); Rednikov 
et. al (1995) assume that a,= ea i, e « 1, and suppose that u depends 
upon a fast variable 9 and a slow time T, such as 



e x = 1, 9t= -V(T), T= et. 
Then equation (1.13) becomes 



0:3 U 000 - V ug + 2ai uug + e [u T + «2 U00 + 0.4 ugg S g + a 5 (u ug) g ] = 0. 

(2.44) 




February 11, 2004 20:10 



WSPC/Book Trim Size for 9in x 6in 



ws-book9x6 



50 



Amplification of Nonlinear Strain Waves in Solids 



The solution of Eq.(2.44) is sought in the form (2.42). In the leading 
order we get 



cxs UQfise — V uo,8 + 2ai uq uq t $ = 0. 



(2.45) 



We are interested in studying localized solutions vanishing together with 
its derivatives at \6\ — » 00 . Then the leading order solution is 

u 0 = — b{T ) 2 cosh ” 2 (b(T)6) (2.46) 

Ot 1 

with V = 4 q !3 b 2 . At order e, we have 



a suites - V Uifi + 2ai (u 0 Ui)g = F, (2.47) 



F = — [uo,t + uo,gg + 0L4 UQfigee + <15 (■ uo uo,e) e ] • 

The operator acting on the function U\ in the lhs of Eq.(2.47) is adjoint to 
that of the Eq.(2.45). Then the secularity condition (2.43) with v t = uq 
yields the equation for b: 



b T = A b 3 (A + Bb 2 ), 



with 



A = « 2 i B = 



4 ( 6a 3 a 5 
7 V «i 




(2.48) 



(2.49) 



The behavior of b (or the sign of br) depends on the signs of A and 
B and on the value of bo = b(T = 0). Indeed, when both A and B 
are positive b diverges while for both negative values it will vanish. For 
A < 0, B > 0 the parameter b vanishes if bo < — A/B while it diverges 

if b 0 > sJ—A/B. The most interesting case occurs when A > 0, B < 0. 
Here b tends to \J— A/B independent of bo- A quantitative description of 
the variation of b can be given. Equation (2.48) may be directly integrated 
giving the implicit dependence of b on T: 



4 B b 2 0 (A + Bb 2 ) 
15 A 2 ln b 2 {A + Bb 2 ) 



4 (b 2 -b 2 ) 
15 Ab 2 b 2 ' 



(2.50) 



We see that b tends to infinity in finite time, T = T* for A > 0, B > 0, 



* 



4 B 



In 



b 2 B 



15A 2 (A + B b 2 ) 



4 

15A &q ’ 



T" 



(2.51) 
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A similar scenario occurs when A < 0, B > 0, bo < \J—A/B. Hence, 
with our asymptotic approximation we can predict when blow-up could 
occur. On the other hand, for both A < 0, B < 0 and A < 0, B > 0, 
the quantity bo > y^—A/B, b vanishes when T tends to infinity. Finally, 
when A > 0, B < 0, the quantity b approaches y/—A/B when T tends 
to infinity, and expression (2.50) provides an analytical description of the 
time-dependent process of the parameter- value selection of the solitary wave 
(2.46), see Fig. 1.17. In the last case, if we additionally assume 03 = 
2040 : 4 / 05 , the asymptotic ’’dissipative” solitary wave (2.46) will tend to 
the exact travelling solitary wave solution ( 1 . 12 ). 

Function u\ ( 6 , T) is a solution of a linear inhomogeneous O.D.E., hence, 
it will contain free parameters depending on T like b(T) in the leading or- 
der problem. Its definition allows to satisfy solvability condition in the 
next order problem and to avoid secular terms in the asymptotic expan- 
sion. Higher order approximations may be studied similarly. The first-order 
solution Mi has been obtained in Rednikov et. al (1995) of the form: 

mi = Aiu 0 ,9 + A 2 6 >cosh ~ 2 (b(T) 9) + A 3 m 0 , 6 ( log(coslr ( 6 (T) 0 ))+ 



Ai i(l — tanh(6(T) 6*)), 

where A, may be found in Rednikov et. al (1995). One can see it predicts 
a plateau behind the solitary wave and does not vanish at minus infinity. 
Uniformly valid solution is found using the matching asymptotic procedure 
may be found in Rednikov et. al (1995). Perturbations of cnoidal wave 
solutions were studied, e.g., in Rednikov et. al (1995); Svendsen and Buhr- 
Hansen (1978). In this case the modulus of the Jacobian elliptic functions 
should be a function of a slow variable, corresponding expressions for the 
derivatives of the elliptic functions with respect to modulus may be found 
in Byrd and Friedman (1954); Newille (1951). 



2.3 Numerical methods 

One can see that both exact nor asymptotic solutions have severe limi- 
tations. In particular they satisfy specific initial conditions and do not 
allow to account for an arbitrary initial disturbance evolution. Therefore 
numerical simulations should be used. Two kinds of nonlinear equations, 
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considered here, may be written either 

u t = L{u) (2.52) 

or 

u tt = L(u) (2.53) 

where L{u) is a certain nonlinear differential operator. Equations of the 
kind (2.52) will be called nonlinear evolution equations while nonlinear 
hyperbolic equations correspond to the class (2.53). Further only nu- 
merical methods are described which are used by the author and his co- 
workers in solving the problems considered in the book. Two main schemes 
are used, pseudo-spectral approach and finite-difference approach. More 
detailed information about numerical modelling of nonlinear wave equa- 
tions may be found, in particular, in Berezin (1987); Dodd et. al (1982); 
Fletcher (1984); Mayer (1995); Sachdev (1987); Zwillenger (1989), see also 
references therein. 



2 . 3.1 Nonlinear evolution equations 

Among the equations of the kind (2.52) main attention has been paid to 
numerical solutions of the KdV equation, see Berezin (1987); Dodd et. al 
(1982) where various difference schemes are discussed. As noted in Berezin 
(1987) implicit three-levels schemes are rather simple and suitable for a 
realization. However, some of the straight finite-difference methods re- 
quire rather small time step for stability, while long-time evolution should 
be studied. It becomes the smaller the higher is the order of a highest 
derivative term in NEE. In particular, it takes place for a simulation of the 
generalized fifth-order KdV equation (1.3), another example has been stud- 
ied in Christov et. al (1997). The problem of an increase of the time step 
may be solved either by a modification of a difference scheme or by using 
a more effective solving procedure, in particular, the fourth-order Runge- 
Kutta finite-difference method. Besides simulation of Eq.(1.3) in Porubov 
et. al (2002), this method was effectively applied to study the evolution 
and interaction of Marangoni-Benard solitary waves governed by Eq.(1.13) 
in Marchant (1996). The same equation has been numerically solved in 
Christov and Velarde (1995) using a four-stage scheme providing third- 
order approximation in time. The implicit predictor-corrector method has 
been applied to nonlinear diffusion equations in Sachdev (1987). 
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Pseudo-spectral methods have been adopted recently, there already exist 
monographs Fletcher (1984); Saclrdev (1987) where they are explained in 
details, see also Zwillenger (1989) and references therein. In a pseudo- 
spectral scheme the space derivatives are approximated very accurately by 
means of the Fourier transforms, in particular, by the fast Fourier transform 
(FFT) algorithm. The fourth-order Runge-Kutta scheme may be used for 
the time derivative Kliaklrandler (1999). It provides very mild stability 
restrictions on the time step. The pseudo-spectral numerical scheme has 
proved very efficient in solving dispersive and dissipative equations. The 
evolution of the initial monotonic solitary wave into radiating or oscillatory 
solitary waves was simulated in Benilov et. al (1993) for the fifth-order 
KdV equation. In a series of papers Salupere et. al (1997); Salupere et. al 
(2001) the solitary wave formation from a periodic input was studied for 
an equation similar to Eq.(1.3), while the paper Salupere et. al (1994) is 
devoted to the KdV soliton detection from a harmonic input. The use of the 
pseudo-spectral approach to account for an evolution of an initial localized 
pulse in framework of Eq.(1.3) is demonstrated in Sec. 1.2.1. Numerical 
solutions of the nonlinear diffusion equations are considered in Fletcher 
(1984); Sachdev (1987). 

Finally, this method allows to study an evidence of the selected solitary 
waves Kliaklrandler et. al (2000) predicted on the basis of the asymptotic 
analysis of Eq.(1.13) in previous section. To reveal the expected behaviour 
of the solutions, the following numerical technique has been employed. 
Since at the very late stage of the evolution the unstable waves are ex- 
pected to be controlled by the finite-amplitude waves found in Kawahara 
(1983), we define the typical wavenumber k c = \J o.i!2a^ corresponding 
to the most unstable linear mode in Eq.(1.13). The corresponding wave- 
length is A c = 2n/k c . The length of spatial domain was chosen to be 256A C , 
i.e. rather long. At the same time, the number of discretization points 
was chosen to be 4096, i.e. A c is covered by 16 points. The latter ensures 
fair resolution of the whole solutions computed. Periodic boundary condi- 
tions have been used for simulations. The pseudo-spectral technique was 
employed for the spatial discretization and the Runge-Kutta fourth order 
scheme for the time advance. The time step was chosen to be 0.01. The 
tests with smaller time steps and better resolution gave indistinguishable 
results. The control of the simulations in the Fourier space shows the good 
resolution of the computed solutions. First, it has been checked that the 
tendency to blow-up, Fig. 2.4(a), the damping, Fig. 2.4(b), and the selection 
have been observed at the a,; values prescribed by the theory in previous 
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(a) Blow-up 




(b) Damping 




distance 



Fig. 2.4 Blow-up (a) and damping (b) of the initial conditions. 

section. In the case of ’’blow-up” , we found that the pulse tends to grow 
rapidly at the time t* rather close to the predicted t* = T* je from (2.51). 
This observation allows to find analytically “time-of-life” of the blowing 
solutions of Eq. (1.13). 

Following Kliaklrandler et. al (2000) let us discuss in detail the selection 
of ’’dissipative” solitary waves occurring at A > 0, B < 0. Choosing the pa- 
rameter values e = 0.1, a\ = 1, = 1, a 3 = 1, <24 = 6/5, as = —2, the 

resulting amplitude of the selected solitary wave is obtained using (2.46), 
(2.49) and b = yJ~A/B. The amplitude, 6013b 2 /a\ = 0.583 and the veloc- 
ity V = 4 a 3 b 2 = 0.389. We consider both the selection occurring from ” be- 
low” when the magnitude of an initial Gaussian pulse is smaller than that of 
the eventually selected solitary waves and the selection from ” above” when 
the selected solitary wave amplitude is smaller than that of the initial pulse 
magnitude. This would permit to separate the selection mechanism from 
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Distance 



Fig. 2.5 ’’Dissipative” solitary wave selection from an initial Gaussian profile with am- 
plitude 0.3 and width 36 units. 



those originating from the growing unstable disturbances Kawahara (1983); 
Oron and Rosenau (1997). One can see in Fig. 2.5 that up to the time 
t ~ 120 an initial Gaussian pulse with the magnitude 0.3 < 0.583 and 
width 36 breaks into a train of three localized pulses aligned in row of 
decreasing magnitude. Due to smallness of e, the influence of the dissipa- 
tive non-KdV terms is small at this stage. This may be seen by compar- 
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Fig. 2.6 ’’Dissipative” solitary wave selection from an initial Gaussian profile with am- 
plitude 1 and width 12 units. 

ison of solutions of Eq.(1.13), shown as solid lines, with pure KdV case, 
«2 = a 4 = a 5 = 0, shown by dashed lines. At later stages the initial pulse 
transforms into a train of the solitary waves. At nonzero e each solitary 
wave amplitude and velocity tend to the values 0.585 and 0.38 in agreement 
with the theory of single solitary wave selection, while each of three KdV 
solitons continues propagation with its own amplitude and velocity. ’’Dis- 
sipative” solitary waves form a bound state Christov and Velarde (1995); 
Nekorkin and Velarde (1994) whose unequal spacing between equally high 
crests reflects the original separation of the solitary waves in the KdV stage 
when higher ideal solitons travel faster. The tail behind the train of soli- 
tary waves appear as a result of short wave instability similar to Kawahara 
(1983). However, the solitary waves have higher velocity than the velocity 
of growing wave packets. As a result, the solitary waves bound state es- 
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capes the destruction induced by the radiation that lags behind. A similar 
finding for the GKS equation was reported in Chang et. al (1995). We see 
on the last stages that the magnitude of the tail saturates. All these struc- 
tures are quite well resolved and similar to those found in Kawahara (1983); 
Oron and Rosenau (1997). They are robust and remain the same under the 
mesh refinement and smaller time steps. 

The selection process realized from ’’above” is shown in Fig. 2.6 when 
initial Gaussian pulse has the magnitude 1 > 0.583 and the width 12. 
Two equal solitary waves with amplitude 0.585 and velocity 0.38 appear 
as a result of decrease of the magnitude of the initial pulse. Again the 
comparison with pure KdV case is shown by dashed lines. All features of 
the selection process are similar to the selection from ’’below”. 

Simulations of Eq.(1.13) with other values of parameters e, Oi such that 
A > 0, B < 0, which are not reported here, show the same features of the 
selection process as described above. 

2 . 3.2 Nonlinear hyperbolic equations 

Among the equations of kind (2.53) we are especially interested in various 
Boussinesq-like long waves equations, that may be written in the form of 
so-called double-dispersive equation Erofeev and Klyueva (2002); Samsonov 
(2001) 

^tt ^1 Uxx ^2 ( U }xx ^3 'Uxxtt T Qq 'U'xxxx = 0- (2.54) 

At q 3 = 0 it corresponds to the classic Boussinesq equation. Besides 
quadratic nonlinear term, <22 (u 2 ) XX i cubic or higlrer-order nonlinearities 
may be considered Christou and Clrristov (2002); Soerensen et. al (1984); 
Soerensen et. al (1987); also higlrer-order derivative terms may be incor- 
porated in Eq.(2.54) Christou and Clrristov (2000); Clrristov et. al (1996). 
A Fourier-Galerkin method were applied to computing localized solutions 
in Christou and Clrristov (2000); Christou and Clrristov (2002). An im- 
plicit difference scheme was developed in Clrristov and Maugin (1995); 
Clrristov and Velarde (1994); Clrristov and Velarde (1995) to account for nu- 
merical solutions of the Boussinesq equation and its generalizations where 
Newton’s quasi-linearization of the nonlinear terms is employed. The so- 
lution in Soerensen et. al (1984); Soerensen et. al (1987) is obtained using 
the Gauss elimination. 

There is an another method developed independently the USSR in 1953 
Godunov and R.yaben’kii (1987); Samarskii and Nikolaev (1989) and in the 
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USA by L.H. Thomas (1949). That is why it is known in the West as the 
Thomas method Morton and Mayers (1994); Richtmyer and Morton (1967). 
As noted in Godunov and Ryaben’kii (1973), this method is justified for 
the solutions of linear problems. However, the author applied it more than 
ten years ago for a numerical solution of Eq.(2.54). One can see further in 
the book that it gives numerical results in a good agreement with analytical 
predictions. Previously the Thomas method was used for both nonlinear 
evolution equationsBerezin (1987) and hyperbolic equationAlexeyev (1999). 
Recently it was successfully applied in Buklranovsky and Samsonov (1998); 
Porubov et. al (1998) for computing rather complicated hyperbolic nonlin- 
ear elastic systems. 

Since the Thomas algorithm is not widely used for solutions of the non- 
linear waves problems, it is useful to describe it in a more elaborate manner. 
Following Godunov and Ryaben’kii (1987); Godunov and Ryaben’kii (1973) 
consider a boundary problem: 

a n u n - 1 + b n u n + c n u n+ 1 = f n , 0 < n < N, (2.55) 

uq = ip, un = ip. (2.56) 

Let us assume 



UO — Li/ 2 Ui + Ki/ 2 , 

where L i/ 2 = 0, Ad/2 = <P- It allows to exclude uq from the equation (2.55) 
at n = 1, giving 



u 0 — L 3 / 2 U 2 + A3/2, 

with L3/2 = —C1/61, Ad/2 = (/1 — cLip)/bi- Following this procedure one 
can obtain 



u n — L n+ i/2«tj+1 + A/i+1/2) 



(2.57) 



with 



r C n fn 1/2 

■ J n+ 1/2 — —J . j > A „+i /2 — t ; , • 

On ~r Un^n- 1/2 0 n + a n Li n - 1/2 



Hence going from 1 to AT one can calculate the coefficients L n+1 / 2 , 
K n +i/ 2 ■ Using (2.57) and the boundary conditions (2.56) at n = N, one 
can obtain un-i- Then coming back from n = N ~ lton = l one can find 
every u n from (2.57) using already known L n+ 1/ 2 , AT n+1 / 2 . It was noted 
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in Godunov and Ryaben’kii (1987); Godunov and Ryaben’kii (1973) that 
small amount of arithmetic operations together with a weak sensitivity to 
the calculation errors are the main advantages of the method. 

Assume u n = u(x n ,tj) is a mesh function of a solution of Eq.(2.54) 
while v n — u(x n ,tj_i), w n = u(x n ,tj_ 2 ) correspond to the mesh functions 
on previous time steps. Then we have for discretization of Eq.(2.54) 

a n — c n — oq, b n Ax 2 ol 3 , 

fn = Ax 2 (w n -2v n )-aiAt 2 (v n -i-2v n + v n+1 )+a 3 (2v n -i-4v n + 2v n+1 - 
w n - 1 + 2w n - w n+ i) - a^At 2 / Ax 2 {v n -2 - 4u n _i + 6v n - 4u„+i+ 
v n + 2 ) - a 2 At 2 {v 2 n _ 1 - 2 v 2 n + v 2 n+1 ), 

where At, Ax are time and space steps respectively. Uniqueness condition, 
\b n \ > \a n \ + \c n | Godunov and Ryaben’kii (1973), is always satisfied at 
positive a 3 . The difference scheme is similar to those used in Soerensen 
et. al (1984); Soerensen et. al (1987) in the case oq = 0 and with higher 
order nonlinearities being taken into account. 



2.4 Use of Mathematica 

Recently various symbolic mathematical programs were developed to pro- 
vide huge analytical operations. One of the most powerful is the Math- 
ematica developed by Wolfram (1999) that is used by the author on all 
stages of his studies. This section is not focused on the detailed description 
of the magnificent abilities of the program. Only some important features 
are considered in relation to the problems studied here. 

Among the advantages of the Mathematica one can mention variety 
of the build-in mathematical functions. In particular, it works efficiently 
with the Weierstrass elliptic function p that may be used for obtaining ex- 
act solutions. Indeed, the procedure described in Sec. 2.1 is based on the 
transformation of the problem of a solution of a PDE to the problem of the 
solution of algebraic equations for the ansatz parameters. These equations 
may be too complicated for manual operations, and use of the Mathemat- 
ica allows to obtain solutions avoiding errors. Design of the Mathematica 
package provides an automatic finding of a solution. An example may be 
found in Parkes and Duffy (1996) where a package is presented to obtain 
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y b 




Fig. 2.7 a) Incorrect and b) correct representations of the function y. 



solutions in terms of the hyperbolic tanh-function. However, it is unlikely 
that an automated method may be applied for finding periodic solutions. 
First, it is necessary to introduce the rules for the elliptic functions deriva- 
tives. Second, sometimes combinations like g 2 C — g% — 4 C 3 or 12C 2 — g 2 
should be kept in the solution for its convenient analysis, see Sec. 2.1.3. 
The Mathematica cannot do it automatically. Finally, it works simulta- 
neously with all possible solutions of the algebraic equations that yields 
huge expressions at the intermediate stages of a solution and may result in 
falling down the evaluation. Hence, it is better to manage the substitution 
procedure step by step manually introducing the commands that provide 
the most efficient line of attack on the problem. 

The algebraic manipulations should be used also for a check of the 
asymptotic solutions and for derivation of the governing nonlinear equa- 




February 11, 2004 20:10 



WSPC/Book Trim Size for 9in x 6in 



ws-book9x6 



Mathematical tools for the governing equations analysis 61 

tions. Numerical abilities may be efficiently used for the analysis of ana- 
lytical relationships, while Fortran, C++ or other languages are preferable 
for numerical simulations of nonlinear PDEs. The Mathematica possesses 
high-level graphic facilities, it admits representation of graphics in many 
formats including PostScript and Encapsulate PostScript. Almost all fig- 
ures in the book are prepared using the Mathematica. Like analytical 
procedures graphics may be automatized. In particular, numerical data 
obtained using Fortran or C-l — I- programs, may be represented by an auto- 
mated procedure, see, e.g., Kliakhandler (1999). 

Sometimes, the Mathematica gives incorrect results. In particular, 
the construction of the profiles shown in Figs. 2.1 and 2.2 requires cal- 
culation of the function y = II \tp, n, k], ip = arcsin(sn(fc x)) expressed 
through the elliptic integral of the third kind. Direct Mathematica com- 
mand Plot [j/ , { 2 ;, 0, 10}] yields the profile shown in Fig. 2.7(a). One can see 
it contains the points where there is no first derivative of y. In order to 
obtain correct smooth profile it is necessary to define y as y = II [ip, n, n], 
where ip is obtained using the Mathematica commands 

ip[x-/] 0 < x < xq\ = ip] ip \x-/\ xo < x < 3 a>o] = 7 t — ip] 
ip[x-/] 3 xq < x < 5 Xq\ = 27 t + ip] 
ip[x-/] 5 Xo < x < 7 £ 0 ] = 37 x — ip] etc. 

where Xq = K(k) / k, K(n) is the complete elliptic integral of the first kind, k 
is the elliptic functions modulus. Then we obtain required smooth profile for 
y shown in Fig. 2.7(b) around which the profiles are developed in Figs. 2.1, 
2 . 2 . 
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Chapter 3 



Strain solitary waves 



in an elastic rod 



The theory of strain waves in solids began to develop over two hundred 
years ago, see about it Love (1927); McNiven and McCoy (1974). During 
the long period only linear theory of elasticity was considered since because 
of the engineering needs and poor experimental facilities. Now the study of 
the material properties Lurie (1990); Murnaghan (1951), acoustic signals 
Biryukov et. al (1991); Oliner (1978); Parker and Maugin (1987) etc. require 
mathematical models based on the nonlinear elasticity. Recent develop- 
ments in general elastic theory may be found in Bland (1960); Lurie (1990); 
Maugin (1993); Maugin (1995); Murnaghan (1951), while nonlinear waves 
in solids were considered in Engelbrecht (1983); Engelbreclrt (1997); 
Erofeev (2002); Erofeev and Klyueva (2002); Jeffrey and Engelbrecht (1994); 
Mayer (1995); Parker (1994); Parker and Maugin (1987); Samsonov (2001). 

In order to go further it is necessary to define a notion of the word 
solid. In macroelasticity solid may be defined as a substance having a defi- 
nite volume and shape and resisting forces that tend to alter its volume or 
shape. From the point of view of the theory of discrete media, solid may 
be considered as a crystalline material in which the constituent atoms are 
arranged in a 3D lattice with certain symmetries. As noted previously 
Jeffrey and Engelbrecht (1994); Samsonov (2001), these definitions comple- 
ment each other allowing to take into account a model of elastic potential 
of atomic interactions, and to cover amorphous, porous or granular media. 

When elastic features are one and the same in any direction we have 
isotropic solids. Sometimes anisotropy is important, e.g., in design of non- 
linear acoustic devices Oliner (1978); Parker and Maugin (1987). Strain 
waves in solids may be classified as follows. Comparing the direction of the 
wave propagation with the particle motion one can distinguish longitudinal 
and shear waves. The former wave propagates along the direction of the 
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particle motion, the latter -perpendicular to it. Waves propagating inside 
solid are called bulk waves. In presence of a lateral surface, surface strain 
waves are possible. In this book main attention is paid to the longitudinal 
bulk waves in isotropic media and wave guides. 



3.1 The sources of nonlinearities 



Among the possible sources of nonlinearity we briefly consider so-called 
geometrical and physical nonlinearities since they affect the strain wave 
propagation to a greater extent. Other kinds of nonlinearities are consid- 
ered in Engelbrecht (1997); Jeffrey and Engelbrecht(1994). The geometrical 
nonlinearity is described by the exact expression of the strain tensor always 
used in the theory of large deformations. Initially the position of a particle 
is accounted for a vector-radius ~P, or an initial or reference configuration 
is defined. Loading forces provide the displacement of particle yielding the 
current or actual configuration characterized by an another vector-radius 
Then the movement is described by the displacement vector, V = 
i? - . However, not only the positions of the particles vary during the 

deformation but also the distances between them. In order to describe the 
alteration of the distance a deformation or strain tensor is introduced. It 
is obtained from the difference between the squares of the arc length in the 
deformed (actual) and undeformed (reference) configuration. In the refer- 
ence configuration the Cauchy-Green finite deformation tensor C is defined 
Lurie (1990), whose general form is 



C = 



V V + (V V) T + v V ■ (V V) T /2 



(written in terms of a vector gradient V V and its transpose (V C) T ). In 
rectangular Cartesian coordinates a;, the components of C may be written 
in a more familiar form 

_ 1 duj du k dui dui 

lk 2 dxk dxi dxi dxk 

where Uk are the components of the displacement vector V. 

As a result of a deformation process there appear stresses. Certainly 
they should be connected with the strains, and this is a source of the phys- 
ical nonlinearity. It was found more than hundred years ago that Hook’s 
linear law of elasticity is insufficient, see about it in Jeffrey and Engel- 
brecht(1994); Samsonov (2001). Since Piola-Kirchoff stress tensor is defined 




February 11, 2004 20:10 



WSPC/Book Trim Size for 9in x 6in 



ws-book9x6 



Strain solitary waves in an elastic rod 65 

in the reference configuration, through the volume density of the internal 
energy II in adiabatic processes (or through the Helmgholtz free energy in 
thermoelastic processes), 

an 

ik ~dCTk 

one can say the physical nonlinearity depends upon the structure of the 
internal (or free) strain energy density n. The energy of deformation must 
be insensitive to the rotation of the reference frame. It was Murnaglran 
(1951) who supposed to develop the energy as a power series in the three 
invariants of the strain tensor, 

n = - 2 nh + - 2 mhh + nh , (3.1) 

where Ik, k = 1,2,3 are the invariants of tensor C: 

Ji(C) = frC, J 2 (C) = [(trC) 2 -frC 2 ]/2, / 3 (C) = detC. (3.2) 

The first two terms in (3.1)account for linear elasticity, hence the second or- 
der elastic moduli, or the Lame coefficients (A, /r), characterize linear elastic 
properties of the isotropic material. Other terms in (3.1) describe material 
or physical nonlinearity Lurie (1990); Murnaglran (1951). Accordingly, the 
third order elastic moduli, or the Murnaglran moduli ( l , m, n) account for 
nonlinear elastic properties of the isotropic material. The energy may be 
written using another set of invariants, 

Ji(C) = trC, J 2 (C) = frC 2 ,J 3 (C) = frC 3 , 

in the form Lurie (1990) 

II = — J 3 + /iJ 2 + — T 3 + rz 2 Ji J 2 + ^3^3- (3-3) 

Since 

Jr(C) =/i(C), J 2 (C) = / 2 (C) — 2/ 2 (C), 

J 3 (C) = J 3 (C) - 3/ 1 (C)/ 2 (C) + 3I 3 (C), 

Eqs.(3.1) and (3.3) coincide if v\ = 21 — 2m + n, 1 S 2 = m — n/2, = n/4. 

In some cases there is a need in more terms in Eq.(3.1): 



TT A “t” 2n 2 T ^ T 2?7i. 3 r T T 

II = — - — Ii — 2 /r/ 2 -f — Ii — 2m/i/ 2 + n/ 3 + 
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Table 3.1 Lame’s and Murnaghan’s modulii, *10 a N/m 2 



material 


A 




l 


m 


n 


Polystyrene 


1.71 


0.95 


-18.9 


-13.3 


-10 


Steel Hecla 37 


111 


82.1 


-459 


-461 


-358 


Aluminium 2S 


57 


27.6 


-299 


-311 


-228 


Pyrex glass 


13.5 


27.5 


14 


92 


420 


S1O2 melted 


15.9 


31.3 


129 


71 


-44 



aJi + CL2I1I2 + (13I1I3 + ( 3 - 4 ) 

The fourth order moduli (ai, ( 12 , 03 , 04 ) can be either positive or negative. 
However, mainly the third-order moduli data may be found in the literature 
Frantsevich et. al (1982); Lurie (1990), some of them are collected in Table 
3.1. 

The series expansions (3.1), (3.4) are convenient to account for the de- 
formation of compressible materials (metals, polymers, etc.), whose yield 
point is small. For incompressible materials, like rubber, there exist an- 
other models, e.g., the Mooney model Lurie (1990). Its generalization, the 
Mooney-Rivlin model, may be used even for compressible materials, its 
application to nonlinear strain waves in a rod may be found, e.g., in Dai 
(1998). 

3.2 Modelling of nonlinear strain waves in a free lateral 
surface elastic rod 

3.2.1 Statement of the problem 

Let us consider an isotropic, axially infinitely extended, free lateral surface 
cylindrical elastic rod, see Fig. 3.1. Axi-symmetry leads to using cylin- 
drical Langrangian coordinates (x,r,ip), where x is the axis of the rod, 
ip e[0,27r], 0 < r < R. When torsions are neglected, the displacement 

vector is V = (u, w, 0). We shall consider the propagation of longitudinal 
strain waves of small but finite amplitude in the rod. Once the reference 
configuration is defined we use Hamilton’s variation principle to obtain the 
governing equations together with the boundary conditions, setting to zero 
the variation of the action functional, 

nt\ nOO nR 

SS = 6 dt2n / dx r Cdr = 0, (3-5) 

J to J —00 J 0 




February 11, 2004 20:10 



WSPC/Book Trim Size for 9in x 6in 



ws-book9x6 



Strain solitary waves in an elastic rod 



67 



r 




Fig. 3.1 Free lateral surface cylindrical rod 



where, C is the Lagrangian density per unit volume, C=K — II, with II 
defined by Eq.(3.1). For the kinetic energy density K we have 



K = 




(3.6) 



where po is the rod material density at t = to- The integration in brackets in 
(3.5) is carried out at the initial time t = to- Initially, the rod is supposed 
to be in its natural, equilibrium state. In absence of torsions non-zero 
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components of the Cauchy-Green deformation tensor C, are 

Cxx ^ x “b f u x ~b VJ X )> C rr w r -b — (u r “b ^C r )> C “h 2^2 ’ 

CV;,; = ^(M r + lKx + U S UI r + W x W r ). 

Hence the invariants in (3.2) may be written as 

-Tl(C) = C xx + C rr + C<p V , 



/ 2 (C) — C xx C rr + C XX C W + C rr C w — C 2 X , 



h(C) = C w {C xx C rr - C 2 rx ). 

The following boundary conditions (b.c.) are imposed for a free lateral 
surface rod: 



w — > 0, 



at r — > 0, 



(3.7) 



P rr = 0 , 



at r = i?, 



(3.8) 



Pvt = 0, 



at r = i?, 



(3.9) 



where the components P rr , P rx of the Piola - Kirclrhoff stress tensor P 
are defined as 

_ / x „ x , w , X + 2n + m 2 

Prr — (A + 2 [i) w r + A — b A u x H — u r + 

r 2 

3A + 6u + 21 + 4m 2 . . u> A + 2/ w 2 

+ _|_ 21) w r I y + 

2 r 2 r z 

, N , x w A + 2/ 9 

(A + 2/) 'iL x vj' l - -\~ (21 — 2m H - Tij Ux — — — — u x 

r 2 



A H- 2/j, + m 9 / \ 

w x + {/i + m) u r w x , 



(3.10) 



_ / , , . x 2A + 2?n — n w 

Prx = H (u r + W x ) + (A + 2/i + m) u r w r H u r — b 

2 r 

. . 2m — n w . . 

( A + 2\i + m) u x u r H w x — b ( /-i + m) w x w r + 

2 r 

(fi + m ) UjiTj . 



(3.11) 
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3 . 2.2 Derivation of the governing equation 

Exception of torsions provides transformation of the initial 3D problem 
into a 2D one. Subsequent simplification is caused by the consideration of 
only long elastic waves with the ratio B/L C 1 between the rod radius R 
and typical wavelength L. The typical elastic strain magnitude B does not 
exceed the yield point of the material. Since Murnaglran’s material have 
small yield points, one can assume B <C 1. The Hamilton principle (3.5) 
yields a set of coupled equations for u and w together with the b.c. (3.8), 
(3.9). To obtain a solution in universal way one usually proceeds to the di- 
mensionless form of the equations and looks for the unknown displacement 
vector components in the form of power series in the small parameters of 
the problem (for example R/ L ), hence, leading to an asymptotic solution 
of the problem. However, this procedure has some disadvantages. In par- 
ticular, comparison of the predictions from the dimensionless solution to 
the experiments suffers from the fact that both B and L, are not well de- 
fined. In particular, solitary wave has an infinite wave length. Further, the 
coefficients of the nonlinear terms usually contain combinations of elastic 
moduli which may be also small in addition to the smallness of B something 
not predicted beforehand. Finally, this procedure gives equations of only 
first order in time, t, while general equations for displacements u and w are 
of the second order in time. Therefore the solution of the model equation 
will not satisfy two independent initial conditions on longitudinal strains 
or displacements Samsonov (2001). 

An alternative is to simplify the problem making some assumptions 
about the behavior of longitudinal and/or shear displacements and/or 
strains in the elastic wave-guide. Referring to the elastic rod these rela- 
tionships give explicit dependence of u and w upon the radius, while their 
variations along the rod axis are described by some unknown function and 
its derivatives along the axis of the rod. Then the application of Hamilton’s 
principle (3.5) yields the governing equation in dimensional form for this 
function. This equation is of the second order of time, hence its solution 
can satisfy two independent initial conditions. Any combinations of elastic 
moduli appear in the coefficients of the equation, hence, subsequent scaling 
may take into account their orders when introducing small parameters. 

For an elastic rod, the simplest assumption is the plane cross section 
hypothesis McNiven and McCoy (1974): the longitudinal deformation pro- 
cess is similar to the beards movement on the thread. Then every cross 
section of the rod remains flat, hence, u = U(x,t) does not change along 
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the radius r. However, this assumption is not enough due to the Poisson 
effect, i.e., longitudinal and shear deformations are related. That is why 
Love (1927) proposed to use a relationship between w and u: w = —r v U x , 
with v the Poisson coefficient. Unfortunately, the plane cross-section hy- 
pothesis and Love’s hypothesis do not satisfy the boundary conditions that 
demand vanishing of both the normal and tangential stresses, P rr and P rx , 
at the lateral surface of the rod with prescribed precision. 

Another theory has been proposed in Porubov and Samsonov (1993) to 
find the relationships between displacement vector components satisfying 
b.c. on the lateral surface of the rod (3.8), (3.9) as well as the condition 
for w (3.7). Later it was developed in Porubov and Velarde (2000). 

Since pure elastic wave are studied, B <C 1, the ’’linear” and ’’non- 
linear” parts of the relationships may be obtained separately. A power 
series approximations is used, as generally done for long wave processes. 
Accordingly, the longitudinal and shear displacement in dimensional form 
are: 

u = ul + unl, ul = u 0 (x,t) + rui(x,t) + r 2 u 2 {x,t) + ..., 
unl = u NL 0 {x,t) +ru NL 1 (x,t) + ..., (3.12) 

w = wl + wnl , u>l = wo(x,t) + rwi(x,t) + r 2 w 2 (x,t) + ...., 
wnl = WNLo(x,t) +rw NL 1 (x,t) + ... (3.13) 

Substituting the linear parts ul and wl (3.12), (3.13) into the b.c. (3.7) and 
in the linear parts of b.c. (3.8), (3.9), and equating to zero terms at equal 
powers of r one obtains uj- and u>k- Using these results the nonlinear parts 
unl, Wnl are similarly obtained from the full b.c. We get Uq(x, t) = U (x, t) 
while other Uk and w k are expressed through U and its derivatives. Then 
(3.12), (3.13) are substituted into (3.1) and (3.6). Running the Hamilton 
principle (3.5) we obtain the single governing equation for the unknown 
function U(x,t). In contrast to the theory based on the plane cross section 
hypothesis and the Love hypothesis, the present theory allows to account 
for nonzero b.c. for the stresses on the lateral surface of the rod. 



3.3 Double-dispersive equation and its solitary wave solu- 
tion 

In order to derive the governing equation for longitudinal strain waves in 
a free lateral surface rod we assume that B ~ R 2 /Li 2 . This assumption 
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provides a balance between nonlinearity and dispersion required for exis- 
tence of the bell-shaped solitary waves of permanent shape. Substituting 
power series (3.12), (3.13) into the boundary conditions and equating to 
zero combinations at each power of the radius, we get 

u = U + a 2 r 2 U xx , (3-14) 

w = h r U x + b 3 r 3 U xxx + B ir U 2 . (3.15) 



where 




2(3 - 2v) ’ 



Bi = 



jz (1 + v) (1 — 2 j /)(1 + v) 
2 + E 



[1(1 — 2v) 2 + 2m(l + v) — nv}. 



Then kinetic and potential energy truncated approximations are 

K = P -^(U 2 + vr 2 [U t U xxt + vU 2 t ] ) , (3.16) 



n = ^ (eU 2 + + vEr 2 U x U xx ^j (3.17) 

Here v and E are the Poisson ratio and the Young modulus correspondingly, 

_ ^ p _ M3 A + 2 /i) 

2(A + fj,) ’ A + n 

while (3 is a nonlinear coefficient, (3 = 3E + 1(1 — 2v) 3 + 4to(1 — 2^)(1 + 
v) + 6 nv 2 . We have to truncate the approximations (3.14), (3.15), (3.16), 
(3.17) in order to be in an agreement with the five-constant Murnaghan 
approximation (3.1), where cubic nonlinear terms are neglected. 

Comparing these relations with those obtained using cross-section and 
Love’s hypothesis, a 2 = 0, b 3 = 0, B 1 = 0, one can see that only the 
term a 2 r 2 U xx makes its contribution into Eqs.(3.16), (3.17), while terms 
b^v 3 U xxx and B 1 r U 2 are needed to satisfy the boundary conditions on the 
lateral surface with prescribed accuracy. 

Substituting (3.16), (3.17) into (3.5) one can obtain the so-called double- 
dispersive equation, DDE, for a strain function v = U x , 

Vtt Qq V x x C^2 ( U } xx Qq V xx tt T Oq Vxxxx ~ 0- (3.18) 
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where 



OL\ 



E f3 

Po 2p 0 



«3 



v(v- 1)R 2 
2 ’ 



Of4 



/ //■ 

2po 



For the first time DDE was derived independently by some authors, in 
framework of the plane cross section and Love’s hypothesis, see about it 
Erofeev and Klyueva (2002); Samsonov (2001). The important difference 
is only in the values of the dispersive terms coefficients «3 and « 4 , in 
Samsonov (1988) it was «3 = ~u 2 R 2 /2, «4 = fiv 2 E 2 /(2p 0 ). 



1 

0.5 




-0.5 

1 

0.5 


iq - * lAaAA ' 150 200 


o o 

CJi I- 1 Cn 


50 100 150 200 


-0.5 

1 

0.5 


50 100 150 200 


-0.5 

1 

0.5 


50 100 150 200 


-0.5 

1 

0.5 


50 100 150 200 


-0.5 

1 

0.5 


50 v 100 v 150 200 


-0.5 

1 

0.5 


50 Tocf 150 200 

0 


-0.5 


50 100 150 200 



Fig. 3.2 Formation of solitary waves from an initial rectangular tensile pulse. 



Eq.(3.18) possesses exact travelling solitary wave solution that may be 
obtained by direct integration, 



i> = Acosh 2 {k(x — Vt)) 



(3.19) 
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with 



A = 



3( Po V 2 -E) )2 



k z = 



E - p 0 V 2 



(3 ’ " 2vR 2 (E + p 0 (v- 1)F 2 ) ' 

Accordingly, V 2 lies either inside the interval 

E 9 9 E 

— < V 2 < c 2 = — , if - 1 < v < 0. 

Po(l ~ v) Po 



(3.20) 



(3.21) 



or in 



cl < V 2 < 77 “ r, if V > 0, 



Po(l - O’ 



(3.22) 




1 

0.5 











- 0.5 


50 


100 


150 


200 



Fig. 3.3 Formation of a train of solitary waves from an initial rectangular tensile pulse. 

When positive Poisson ratio is positive, more precisely, 1/2 > v > 0, 
the sign of the amplitude is defined by the sign of f3. The nonlinearity coef- 
ficient is the only coefficient carrying an information about the Murnaghan 
moduli. Hence, they define whether tensile or compression solitary wave 
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may propagate for a given elastic material of the rod. We see that it is the 
mixed dispersive term a 3 v xxtt , who establishes the permitted finite interval 
for the wave velocity. 

Exact solution requires specific initial conditions. It is of interest to 
know how rather arbitrary localized initial pulse evolves. To this purpose 
numerical simulation of the equation has been performed using numerical 
method explained in Sec. 2.3.2. 
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Fig. 3.4 Delocalization of an initial compression pulse. 



It was found that rather arbitrary initial pulse splits into the train of 
solitary waves or evolves into small amplitude oscillating wave-packet ac- 
cording to the predictions about wave parameters done on the basis of exact 
single travelling solitary wave solution. Indeed, assume the material elastic 
features of the rod yield (3 > 0. Shown in Fig. 3.2 is the evolution of an 
initial tensile rectangular pulse into single solitary wave, while massive ini- 
tial pulse splits into the sequence of solitary waves, see Fig. 3.3. However, 
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each solitary wave is accounted for the exact single travelling wave solution 
(3.19). Figure 3.4 demonstrates no solitary wave generation from the initial 
compression pulse. However, when the material of the rod possesses 0 < 0, 
only compressive solitary wave arises from an initial compression localized 
input. To sum up, single travelling wave exact solution provides correct 
predictions about solitary waves formation in an unsteady process of an 
arbitrary localized input evolution. 

Finally, let us estimate the deviation from the plane cross section caused 
by the second terms in the expression (3.14). Indeed, the curvature h 
of the profile of longitudinal strain u along the radius of the rod is h = 
| i( rr | /(I + n 2 ) 3 / 2 = \vU xx \ /(I + jz 2 r 2 /7 2 a ,) 3 / 2 . It was found in Porubov and 
Samsonov (1993)that the variation of longitudinal strain u along the radius 
is negligibly small when U x is the solitary wave (3.19) that is important for 
the experimental observation of the waves. 



3.4 Observation of longitudinal strain solitary waves 

We briefly consider the recent successful experiments on solitary wave ob- 
servation in a transparent rod, details may be found in Dreiden et. al (1995); 
Samsonov (2001). The optical methods were used because they are prefer- 
able to study transparent optical phase inhomogeneities. They allow not 
only to visualize inhomogeneity but also to determine its parameters, and 
on the other side, being contactless, they do not introduce any disturbances 
in an object under study. All optical methods record the alterations of re- 
fractive index in an object, when studying optically transparent phase in- 
homogeneities. Shadowgraphy is more convenient to record a considerable 
refractive index gradient, for example, caused by a strong shock wave prop- 
agation. It was shown theoretically in our case that a strain solitary wave 
is a propagating long density wave of small amplitude. Interferometry is 
the most appropriate for such waves study because it allows to observe and 
measure with sufficient accuracy even small refractive indices alterations. 

Holographic interferometry method, used in the experiments Dreiden 
et. al (1995), has several advantages in comparison with the conventional 
optical interferometry. In particular, limitations to the optics quality are 
considerably lower because wave fronts to be compared pass through the 
same optical path. For this reason both waves are distorted to the same 
extent and possible defects in optical elements and experimental cell do not 
affect the resulting interference pattern. However, the choice of an optical 
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recording method allows to study, in general, only elastic materials, that 
are transparent for the given light wave length. 

To check that the excited strain wave possesses indeed the solitary wave 
feature to conserve its shape, it is necessary to follow in observations a 
propagation along an extended elastic wave guide. However, in a wave 
guide made of material highly absorbing linear elastic waves, it will be 
sufficient to detect the constant shape wave propagation at much shorter 
distance. 

Based on the results of the analysis presented above the transparent 
polystyrene SD-3 has been chosen as an appropriate material for a wave 
guide manufacturing. The elastic properties of it are given by a set of 
parameters v = 0.35; (3 = — 6 • 10 10 N/m 2 ; c* = 1.8 • 10 3 m/sec, see 

Frantsevich et. al (1982). 

The experimental set-up used to generate and observe the strain soli- 
tons, consists of a basin where the rod is submerged into the water, a 
device to produce the initial shock wave, a holographic interferometer for 
the recording of a wave pattern, a synchronizer and a laser radiation energy 
meter. The waves inside the rod are generated from a primary shock wave 
produced in the water near the edge of the rod by laser vaporization of a 
metal target. The first exposure of the hologram is carried out to obtain 
the hologram of undisturbed wave guide. The second exposure is produced 
by a laser pulse synchronized with the prescribed stage of the wave propa- 
gation. Observations are carried out in the transversal direction, and two 
cut-off were made parallel to the rod axis in order to make transparent 
the central part of the rod. The carrier fringes on interferograms, obtained 
due to the reconstruction of doubly exposed holograms, occur due to the 
wedge turn between the exposures. The longitudinal strain wave patterns 
are recorded at various distances from the input edge of the rod, that are 
attained by the cell displacement along the axis of wave propagation. 

The solitary wave parameters are calculated based on the data of the 
holography interferograms obtained. Typical interferogram is shown in 
Fig .3.5. Only central part of the PS rod is transparent thanks to the 
vertical cuts. Curved interference fringe is extracted from the area inside 
the rod and placed below the interferogram. Note that the interferometric 
pattern does not exhibit a standard bell-shaped image of a shallow water 
solitary wave since the strain solitary wave is, in fact, a longitudinal density 
wave in solids. 

The solitary wave amplitude can be calculated using the interferometry 
fringe shift A K measured in the interference pattern, see Fig. 3.5. Let 2 h 
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Fig. 3.5 Experimental observation of longitudinal strain solitary wave in an elastic 
rod(after Dreiden et. al (1995)). 



be a distance passed by the recording light across the rod, i.e. , precisely, 
the distance between two longitudinal cut-off. Before the deformation the 
phase variation A<pi of the light wave having the length A, is caused by the 
laser light propagation along the distance q — 2h through the water and 
the distance 2 h through the rod ( where q is the distance between the cell 
walls) Samsonov et. al (1998): 



no(q — 2h) + 2hni = — A<f>i. (3.23) 

27T 

Here no,ni are the refraction indices of water and the elastic material 
before deformation. After the deformation the refraction index value of 
the rod changes to « 2 - Moreover, the distance, which light passes through 
the rod and water, varies due to the deformation of the rod. As a result 
we obtain the formula for the magnitude of the light wave phase variation 

A(f>2- 



A A , 

= 



no(q — 2h — 2 Ah) + 712(2 h + 2 Ah) 



(3.24) 
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Evidently, the interferometry fringe shift A K is defined as 



A K 



A0 2 — A0i 

27 r 



(3.25) 



The new value of the refraction index of the deformed rod ?i 2 is caused by 
the local density variation: 



^ (3.26) 

p ni - 1 

On the other side, one can obtain the density variation from the solution 
of a static linear problem on uniaxial compression (or tension), see Lurie 
(1990), and as a result the following relationship is valid: 

— = U x (2u — 1). (3.27) 

P 

Then we get finally from (3.25): 



U x 



XAK 

2h[{n\ — 1)(1 — 2v) + v{ni — ?ro)] 



(3.28) 



The amplitude is determined by the maximal fringe shift value. Deriva- 
tion of the relationship (3.28) shows that the length L of the solitary pulse 
may be directly determined from the interferogram as the length of the 
fringe shift perturbation between two undisturbed areas. 

In our experiments Dreiden et. al (1995) we used a long polystyrene rod, 
149 mm long. The nonlinearity parameter (3 < 0 for polystyrene, and we 
can anticipate the appearance of only compression solitary wave. It is fairly 
long and keeps its shape on propagation. As follows from (3.28) the wave 
shown under interferogram in Fig. 3.5 is indeed a longitudinal compression 
wave. There is no any tensile wave around the observed wave, moreover, in 
the water surrounding the rod the interference patterns remains horizontal, 
i.e., undisturbed, that proves the wave that is measured propagates inside 
the rod. Therefore, the theory developed for an infinite cylindrical rod 
allows to predict an existence of solitary waves in a finite-length rod with 
vertical cuts. 

Unfortunately only the wave amplitude may be measured more or less 
precisely, hence full quantitative comparison with the theory is impossible. 
The unsteady process of the solitary wave generation cannot be accounted 
for the theory. The propagation of a shock wave before the solitary wave 
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contradicts numerical results shown in Fig. 3.2. For a long time solitary 
waves were observed only in a wave guide made of polystyrene. Very re- 
cently they were generated also in a plexiglas rod Samsonov et. al (2003) . 



3.5 Reflection of solitary wave from the edge of the rod 



Following Dreiden et. al (2001) let us consider a semi-infinite homogeneous 
rod, —oo < x < X. The procedure developed in section 3.2 may be applied 
but now the elementary work done by the external forces should be included 
into the Hamilton principle, and the integration with respect to x is carried 
out over — oo < x < X, 



rti 


j-X fR 


/ dt 


2n / dx r Cdr 


't 0 


J — oo J 0 



SAdt = 0. 



'to 



(3.29) 



The Lagrangian density per unit volume, C=K — n, is defined as before. 
According to Lurie (1990) we cannot call Eq.(3.29) variation principle since 
external forces are not necessary potential ones. Use of the Hamilton prin- 
ciple yields the DDE (3.18) as a governing equation for longitudinal waves 
together with the boundary conditions requiring zero values for v and its 
derivatives as x — > — oo, while on the rod end, x = X, they depend upon 
the type of clamping. If the rod end is assumed free, the elementary work 
of external forces at the end of the rod, x = X, is zero, SA = 0, and Eq. 
(3.29) yields 



v = 0, v xx = 0. 



(3.30) 



When the end is clamped, the elementary work is not determined, but 
obvious kinematic reasons require zero displacement and its velocity, 

U = 0,u t = 0 , (3.31) 

that may be rewritten in terms of strains, 

v x = 0, v xt = 0. (3.32) 



Since DDE does not possess an exact solution describing interacting waves 
moving in opposite directions, it is necessary to apply an asymptotic tech- 
nique, reducing DDE to a nondimensional form. Introducing the scales, L 
for x, L/co for t, and B/L for v, where cq = \/E/ p 0 is the so-called rod 
velocity. A small parameter e is chosen as e = B = Ft? /L 2 « 1, to 
balance nonlinearity and dispersion. 
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80 Amplification of Nonlinear Strain Waves in Solids 

Assume the asymptotic solution v depends, in addition to x, t, on a 
slow time r = et. The solution is sought in power series in e: 

v = vq + ev i + . . . . (3.33) 

Substituting (3.33) into dimensionless DDE, and equating the terms with 
the same power of e we obtain the D’Alembert solution for Vq: 

v 0 = v 0 i (0, t) + v 02 (3.34) 

where 9 = x + t, ip = x — t. At order e there is a linear equation for V\, 



2vi,8iI> = 2l>oi,0 T — 2 uo2,v>t 



2 E 









-^-(voi,8eoe + fo2,-0V’bb) (3.35) 

The absence of secular terms leads to two uncoupled KdV equations for the 
functions vq\ and U 02 , 



2voi,t — ^(uoi )o — 2~vqi,889 = 0, (3.36) 

P v 2 

2uq2,t T ^j^{Vo2)ip 4 ^1 , 02 ,tptptp ~ (3.37) 



and the bounded solution of Eq. (3.35) is 

vi = — U 01 U 02 + vu(9, t) + r). (3.38) 

h/ 

Substituting the soliton solutions of the KdV equations (see Chapter 1) 
into the leading order solution (3.34), we obtain 

6^2 

Vq = — - — k 2 [cosh -2 (k[x + (1 + ev 2 k 2 )t — xoi) + 

P 

cosh -2 ( k[x — (1 + £v 2 k 2 )t — xo 2 ])]- (3.39) 

where X 01 are the constant phase shifts. It follows from (3.39) that the type 
of the strain wave depends upon the sign of /3 like in previous section. Both 
conditions (3.32) on the clamped end are satisfied if 202 = 2X — xqi- In this 
case reflection of the solitary wave occurs with no change of its shape. On 
the contrary, one cannot satisfy boundary conditions (3.30) at a free end of 
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Fig. 3.6 Reflection of the solitary wave from a free end of the rod. 



the rod. It means that the reflected solitary wave does not propagate but 
disappears due to dispersion. 

It may be clearly seen from numerical simulation of the reflection Drei- 
den et. al (2001). Again the implicit finite-difference scheme explained in 
Sec. 2.3.2, is used. Implementation of the boundary conditions at the end 
of the rod is effected by means of symmetric continuation of the calculation 
area beyond the real rod end, so that the area of calculation would occupy 
the interval 0 < x < 2X. In case of free end of the rod the sign of /? in 
the interval X < x < 2X is chosen opposite to that used in the interval 
0 < x < X. Initial pulses are assumed to be the equal magnitude KdV 
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Fig. 3.7 Reflection of the solitary wave from a fixed end of the rod. 

solitons located centrally symmetric relative to the genuine end x = X. 
The type of an initial solitary wave depends upon the sign of (3. In par- 
ticular, if (3 > 0 the initial solitary wave in the interval 0 < x < X is a 
tensile wave, while in the interval 0 < x < 2X we impose the same tensile 
solitary wave for the clamped end or the compression solitary wave for the 
free end. Initial velocities of the solitary waves are chosen equal and taken 
to be directed towards each other. 

The reflection from the free end is shown in Fig. 3.6. The right-hand 
side of the figure corresponds to the free end. It is seen that the amplitude 
of the solitary wave propagating from left to right decreases as it reaches 
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Fig. 3.8 Observation of a reflection of the solitary wave from a free end of the rod. 



the end of the rod. The reflected wave is of opposite, first, its amplitude 
grows, then the wave is dispersed, and no localized strain wave is observed 
near the input end of the rod. 

In case of the clamped end numerical results are shown in Fig. 3.7. In 
agreement with the theory the amplitude of the solitary wave is seen to 
be nearly twice as large when it reaches the end of the rod. The reflected 
wave has the same amplitude and velocity as the incident one. Keeping its 
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Fig. 3.9 Reflection of the solitary wave from a clamped end of the rod. 



shape, the reflected wave propagates towards the input end of the rod. 

Experimental observation of the solitary strain wave reflection were car- 
ried out in Dreiden et. al (2001) using the same technique as for the study 
of the solitary wave propagation. Shown in Fig. 3.8 is the reflection from 
the free end. Footnotes in Fig. 3.8(a) and Fig. 3.8(b) demonstrate the 
decrease of the amplitude of the initially generated compression localized 
wave shown in Fig. 3.8(a), cf. with the first and the fifth stages in Fig. 3.7. 
It was found in experiments the absence of any reflected localized waves 
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Fig. 3.10 Reflected wave comes back to the input end of the rod. 



moving towards the input end of the rod, cf. with the last stage of Fig. 3.6. 
Hence all observations verify theoretical predictions. 

Figures 3.9, 3.10 demonstrate the reflection of the strain solitary wave 
from the rod end attached to the right to a massive brass plate. The 
material of the plate is chosen so as to avoid penetration of the wave energy 
outside the rod. It was found experimentally in Dreiden et. al (2001) that 
the acoustic resistance of the plate should be much greater than that of 
the material of the rod (polystyrene). It is seen that the amplitude of the 
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incident solitary wave, Fig. 3.9(a), is almost doubled by the fixed end of 
the rod, Fig. 3.9(b)in an agreement with the first and the fifth stages in 
Fig. 3.7. Figure 3.10 shows the reflected solitary wave moving to the left 
at a distance of 140 mm from the fixed end. It has the same amplitude 
and velocity as the incident one, cf. the last stages in Fig. 3.7. In order to 
see it is necessary to subtract the fringe shift outside the rod (marked by 
B in Fig. 3.10) from the fringe shift measured inside the rod (marked by 
A in Fig. 3.10). The wave observed in experiments demonstrates the main 
feature of the solitary waves to keep their shapes after collisions, hence it 
confirms also the fact that the observed localized wave is indeed the strain 
solitary wave predicted by the theory. 
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Chapter 4 

Amplification of strain waves in 
absence of external energy influx 



As we seen in the previous chapter, strain solitary wave propagates without 
change of the shape in an uniform rod with a free lateral surface. However, 
the wave may still exist even in presence of an external medium or when 
a material of the rod is microstructured. At the same time, the shape of 
the wave varies when the balance between nonlinearity and dispersion is 
destroyed. The simplest reason is the varying cross section of the rod. 



4.1 Longitudinal strain solitary wave amplification in a nar- 
rowing elastic rod 

The section is devoted to the theoretical and experimental description of 
the propagation and amplification of the strain solitary wave (soliton) in 
a cylindrical nonlinearly elastic rod with varying cross section. We call it 
an inhomogeneous rod in the following for convenience, while the rod with 
permanent cross section will be called the homogeneous one. We follow 
Samsonov et. al (1998) where the results presented below were published 
for the first time. 



4.1.1 Governing equation for longitudinal strain waves 
propagation 

Let us consider the wave propagation problem for an isotropic infinite non- 
linearly elastic compressible rod with varying cross section, see Fig. 4.1. 
Introducing the cylindrical Lagrangian coordinate system (x, r, ip), where x 
is the axis along the rod, ip e[0, 2 tt] , 0 < r < R( x) < Rq, Rq— const, one 
can write the displacement vector V = ( u , w , 0), if torsions are neglected. 
Basic equations, describing the nonlinear wave propagation in the initial 
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Fig. 4.1 Free lateral surface cylindrical rod with varying cross section. 

configuration, are obtained from the Hamilton variation principle (3.5). 

We are studying long nonlinear longitudinal strain waves (density 
waves), that allows to do some simplifications, namely, the relationships 
between longitudinal and transversal displacements u and w. To find them 
one needs to satisfy the boundary conditions on the free lateral rod surface 
r = R(x) , namely, the absence of both the normal and tangential stresses 
at every moment. We introduce the small parameter e, taking into account 
that the waves under study should be elastic waves with sufficiently small 
magnitude B, B « 1, as well as sufficiently long waves with the length 
L , so as the ratio Rq/L « 1, where Rq is the maximal value of r(x) 
along the rod . The most important case occurs when both nonlinear and 
dispersive features are in balance and small enough: 

£ = B = ( x) « L ( 41 ) 

The boundary conditions have the form (3.7)-(3.9). The unknown func- 
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tions u, w are expanded in a power series (3.12), (3.13). Substituting 
them into the boundary conditions and following the procedure explained 
in Sec. 3.2, one can obtain 



u L = U(x,t) +r 2 -U xx , w L = -r vU x - 



2(3 - 2v) 



U x 



(4.2) 



WNL = r 



v{l + v) {l~2v){l + v) 
2 + E 



(1(1 - 2u) 2 + 2m(l + v) — nv) 



U', 

J (4.3) 

Higher order terms may be obtained by the same way, however, they will 
be omitted here because of no influence on the final model equation for 
the longitudinal strain waves equation to be derived using the Hamilton 
principle (3.5). Substituting expansions (3.12), (3.13) with (4.2), (4.3) into 
the formulas for K and n, one can find respectively: 



K = f ([/? 



+ vr 2 [U t JJ xx t + ) , 



(4.4) 



n = ; I EUl + ^ul + uEr 2 U x Uxxx 



(4.5) 



where /? = 3E + 2l(l — 2j /) 3 + 4to(1 + ^) 2 (1 — 2^) + 6?w 2 becomes the only 
coefficient depending on the nonlinear elasticity of the rod, it coincides 
with that obtained in Sec. 3.2.3. Substituting (4.4), (4.5) into (3.5) and 
calculating SS = 0, one can obtain the following nonlinear equation: 



U tt - 



^ d_ 

R 2 dx 



[R 2 U X = 



d ' 

R 2 dx 



P 



v d 



^ r(, '- 4<S (ru “ ) + T r(/ - 



R 2 dx 



4 



R a U t 



d?_ 

dx 2 



(R 4 U x 



vR 2 



U T 



(4.6) 



where c* is the so-called ’’rod” wave velocity, c 2 = E/p 0 - 

Let us consider now the rod which cross section varies slowly along the 
a:— axis, which is described by a function R = R("/x), 7 << 1. Introducing 
the notation: v = U x , r = tc* and differentiating the equation (4.6) on 
x, we obtain an equation 



did 



Vtt dx R 2 dx 



R 2 V + (^) + aR^Vrr - bR^Vxx - 4 bR i R x v x J = 0 

(4-7) 
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with the dispersion terms coefficients a and b , a = — [i/(l — i/)]/2, b = 
—v 1 2. Two first terms here describe a common linear wave, the third 
governs the nonlinearity, and the two following terms are responsible for 
dispersive features of the rod, while the last term, being of the same order, 
looks like a dissipative one, it occurs due to the cross section variation. 



4 . 1.2 Evolution of asymmetric strain solitary wave 

To describe the evolution of a travelling strain wave v we introduce the 
phase variable 9 and the slow variable X = jx, as follows: 

9r = -1, 9 X = A(X). (4.8) 

Using the asymptotic method explained in Sec. 2.2, the solution of Eq. 
(4.7) is obtained in new variables in the power series in 7: 

v = v 0 + 7U1 + (4.9) 

Substitution (4.9) into the equation (4.7) gives in leading order of 7 the 
ODE reduction of the nonlinear double dispersive equation (3.18) for i>o: 

(1 ^ A 2 )vQfiQ — A 2 0)00 — R 2 A 4 (a — b)v o,eeee = 0, 
whose solitary wave solution is Samsonov (2001): 



v 0 = ~j a cosh" 2 (k(X) [9 -6 0 (X)]), (4.10) 



depending upon the varying parameter a = a(X), a > 0, while A, k are 
expressed through it: 






1 2 a (l + a ) 

1 + q’ ' 4i? 2 [a(l + a) — b\ 



(4.11) 



Since the parameters of the solitary wave depends upon the slow coordinate, 
the wave is asymmetric like shown in Fig. 1.18. 

Both A and k will be real in Eq.(4.11) for most standard elastic materials 
(having the Poisson coefficient v > 0) if the value of the function a is inside 
an interval : 



v 



0 < a < 



1 - v' 



(4.12) 
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Then the type of the strain wave (4.10) (compression or tensile one) is 
defined only by the sign of the nonlinear coefficient P, which depends on 
the elasticity of the rod material, respectively. 

Let us study a distortion of the solitary strain wave due to the ” geomet- 
rical” inhomogeneity considered. Indeed an inmogeneous linear equation 
holds for Vi at order 0(j) 



(1 - A 2 )vi,gg - A 2 j^(viV 0 )gg - R 2 A 4 (a - b)viMggg = F, (4.13) 

where 

F = —j^~A{ 2uo,0 + + vR 2 [3(l> — 1) + 5A 2 ]vo^ggg) + 

Ax(vo t g + ^^(r'o)s + vR 2 [ — — 1- 3A~\vQjgg) + 

A(2v 0i gx + ~^{vo)gx + vR 2 [ — — 1 - 2 A 2 ]vQfiggx)- 

Then the equation for the amplitude a arises from the secular term absence 
condition (2.43): 



In 



R 2 a 2 \ 4 bk 2 R 2 A 4 

2 kA 3 ) „ + 5 



(In 2 R 4 a 2 Ak) x = 0, 



(4.14) 



that after use of Eq.(4.11) is reduced to a nonlinear first order ODE for an 
amplitude variation 



Rx _ ( 1 1 1 1 

Al ~ ax \6{l- D + a ) "" 2a _ 3(1 -D 1 +a) _ 3(1 - D 2 + a) 

(4.15) 

where: 

b _ i _ 2±yw^ 

a (1 — v) ’ 1,2 5(1 — v) 

Taking the restrictions for a, Eq.(4.12), into account, we conclude that the 
expansion in the brackets on the right- hand side of Eq.(4.15) is always 
positive. Therefore the unbounded growth of the amplitude of the solitary 
wave occurs with the radius decrease, while the increase of the radius is 
accompanied by the decay of the amplitude. When the radius no longer 
alters at X > X*, see Fig. 4.1, the amplitude and the velocity of the solitary 
wave remain constant. The solitary wave evolution in this case is similar 
to that shown in Fig. 1.19. However, there is no tendency to a finite value 




